﻿Capitolul Ecuații diferențiale de ordinul întâi rezolvabile prin metode elementare Definiția O ecuație diferențială de ordinul întâi este o relație de dependență funcțională de forma g(t, x, x) = ( ) înțre funcția idențică t t definița pe ințervalul I necunoscuț, o funcție necunoscuță x si derivața ei X definițe pe același ințerval în ecuația ( ) funcția g se considera cunoscuta, iar rezolvarea ecuației înseamnă determinarea funcțiilor necunoscute x care verifica ecuația Definitia O funcție reala x de clasa C definiță pe un ințerval deschis I C IR se numesțe soluție a ecuației ( ) daca pențru orice t G I țriplețul (t,x(t),x(t)), aparține domeniului de definiție al lui g si g(t,x(t),x(t)) = ( ) Graficul unei soluții: Г = {(t,x(t))|t G I} se numesțe curba ințegrala La începuț vom prezența cațeva cazuri particulare de asemenea ecuații, care se rezolva cu mețode elemențare si probleme concrețe din diferițe domenii care au condus la asemenea ecuații CAPITOLUL Problema primitivei Ecuații diferențiale de forma x = f( Problema O conducta termica are diametrul (cm) și este izolata cu un strat cilindric de (cm) grosime Temperatura conductei este (°C), iar temperatura mediului exterior este (°C) i) Care este legea de variatie a temperaturii in stratul izolant în cazul staționar? ii) Ce cantitate de caldura cedează fiecare metru de conducta în ore? Se da coeficientul de conductivitate termica k = , (W/K • m) Rezolvare: Fie t distanta unui punct din stratul izolant si axa conductei termice, t G ( , ) x - m si x(t) temperatura în acest punct Temperatura este functia necunoscuta si depinde de t, iar functia x = x(t) descrie variatia temperaturiii în stratul izolant i) Pentru determinarea functiei x(t) folosim legea lui Fourier: cantitatea de caldura Q cedata în unitatea de timp în regim stationar pe suprafata laterala a cilindrului de raza t este proportionala cu produsul dintre aria laterala a cilindrului si variatia temperaturii dx: -k's =Q ( ) unde k este conductivitatea termicăa a materialului izolant Aria laterala a cilindrului de rază t si lungime l, este S(t) = n • t • l Rezulta: dx Q dt n kl t ( ) Prin urmare avem de determinat o functie x(t) care veifica ( ) Din teoria primitivelor rezulta ca orice functie x(t) care verifica ( ) este datăa de formula x(t) = - ănRkz'ln г + C ( - ) în care t G ( , ) si C este o constanta oarecare Determinarea legii de variatie cerute revine la selectionarea acelei funcții x(t) din familia ( ) care Problema primitivei Ecuații diferențiale de forma x = f (t) verifica condițiile: pentru ti = (cm) avem x(t ) = (°C), și pentru t = (cm), x(t ) = (°C) Impunând aceste condiții, rezulta ln o- C = + -= , (°K) si ln Q n k l In , de unde x(t) = , - , In t ii) Folosind valorile numerice rezulta Q iar pentru cantitatea de caldura cedata de fiecare metru liniar (l= m) in ore, q = ® • • (s) Trecem acum la cazul general de rezolvare a unei ecuatii diferentiale de forma X = f (t) in care f este o functie reala continua definită pe un interval I C IR , considerata cunoscuta Din teoria primitivelor se §tie ca, daca f este o functie reala continua definita pe un interval I C IR , atunci exista o familie de functii reale de clasa C definite pe I a caror derivata este functia f Aceste functii difera între ele printr-o constanta si se obtin cu formula: x(t) = /* f (t)dr + C ( ) Jt* în care C este o constanta reala iar / f (t)dT este o primitiva a functiei f t* Pentru t € I si xo G IR exista o singura solutie x = x(t) a ecuatiei ( ) care verifica conditia x(t ) = x si aceasta este data de formula x(t) = x + / f(s)ds Jt ( ) Problema determinarii acelei solutii a ecuatiei x = f (t) care verifica conditia x(t ) = x se numeste problema cu date initiale sau problema Cauchy Intr-o asemenea problema t si x se considera cunoscute si se numesc date initiale Problema în sine se noteaza traditional astfel: x = f (t) x(t ) = x ( ) si solutia ei cu x(t; t , x ) CAPITOLUL Concluzii Exista probleme de fizica care conduc la ecuații diferențiale de forma x = f (t) (numita problema primitivei) în care f este o funcție reala continua definita pe un interval deschis (a,b) C IR Oricare ar fi solutia x = x(t) a ecuatiei diferentiale x = f (t) exista o constanta reală C astfel încat x(t) = [ f (t)dr + C, (V)t G (a,b) Jt* Oricare ar fi t G (a, b) si x G IR exista o singura functie x = x(t) definita pe (a, b) care este solutia problemei cu date iiihiale x = f (t) x(to) = xo Exerciții Sa se determine solutiile urmatoarelor ecuatii diferentiale (cu calculatorul): t t a) x = + t + t ; t G IR R : x(t) = — + — + t + C b) x = -; t > R: x(t) = lnt + C c) x = + sin t + cos t; t G IR R: x(t) = t — cos t + sin t + C d) x = ; t G R>l R: x(t) = arctan t + C — t e) x = ;t G (— ) R: x(t) = lni+t+ C f) x t ; t G IR [— , ] R: x(t) = ln(t + ) + C Problema primitivei Ecuații diferențiale de forma x = f (t) g) x = e t + sin t; t G IR R: x(t) = e t — cos t + C h) x = et ; t G IR Sa se rezolve următoarele soluțiile (cu calculatorul): R: se determina numeric o primitiva a lui et , de exemplu / es ds probleme Cauchy si sa se reprezinte grafic a) x = + t + t , t G IR , x( ) = t t R: x(t) = y + +t+ b) x = -, t > , t x( ) = R: x(t) = ln t c) x= + sint + cos t, tGIR , x(—n) = R: x(t) = —cos t + ^sin t+ + +n d) x , t G IR , + x(- ) = - R: x(t) = arctan t + - n — - Ț •« , t— t R: x(()=lnv • ■ — '' x( ) = R: x(t) =l^ + ^/ + j+ln + CAPITOLUL Ecuații diferențiale autonome x = g(x) Problema O racheta meteorologica este lansata vertical în sus cu viteza inițiala de (m/s) Rezistenta aerului frâneaza miscarea ei și-i comunica acceleratia — k • v (t), v(t) fiind viteza rachetei la momentul t iar k o constanta pozitiva i) Sa se afle timpul in care racheta atinge inaltimea maxima ii) Sa se afle înalțimea maxima la care se ridica racheta Rezolvare: i) Acceleratia totala a rachetei, in lansarea pe verticala în sus este a = — (g + kv ) unde g ~ (m/s ) este acceleratia gravitațională, iar k o constanta pozitiva considerata cunoscuta Legea de miscare a rachetei se scrie astfel: dv dt (g + kv ) ( ) Funcția v care intervine în ( ) reprezinta viteza rachetei si este necunoscuta Ea trebuie gasita pentru ca apoi egaland-o cu zero (aceasta înseamna ca racheta a atins înălțimea maxima) sa găsim timpul în care racheta atinge înalțimea maxima Din ( ) si din inegalitatea g + kv > rezulta egalitatea g + kv dv dt Trecand la primitive se obține egalitatea g + k v (t ) dv —dr = dr din care printr-o schimbare de variabilăa rezultăa — kT\t- sau v(t) = • tanf( — kt + C) Ecuații diferențiale autonome x = g(x) Constanta C se determina din condiția inițială v( ) = (m/s) și se obține iar timpul ti dupa care racheta ajunge la înălțimea maxima se determina din condiția v(t ) = si se obtine ti = (s) ii) Pentru a gasi înalțimea maxima la care se ridica racheta se noteaza cu x(t) înalțimea la care se afla racheta la momentul t Funcția x(t) este necunoscuta si pentru determinarea ei se ține seama ca viteza v(t) este derivata funcției x(t): si ca x( ) = (racheta pleaca de pe sol) Determinarea funcției x(t) care verifica aceste condiții este o problemă Cauchy de forma x = f (t), x(t ) = x si rezolvarea ei a fost facuta în § Se determina soluția x(t; , x ) a problemei Cauchy si se calculeaza apoi x(t ; , x ) Aceasta este înalțimea maxima la care se ridicăa racheta meteorologicăa Raționamentul prezentat la rezolvarea punctului i) al problemei poate fi generalizat pentru determinarea soluțiilor unei ecuații diferențiale de forma: x = g(x) ( ) în care g este o funcție reala continua definita pe un interval J C IR , care nu se anuleaza (g(x) = (V)x G J) si este cunoscuta Intr-adevar, daca o funcție reala x : I J este o soluție a ecuației ( ) atunci pentru orice t G I avem dx dt g(x(t)) sau g(x(t)) dx dt CAPITOLUL Trecând la primitive rezulta egalitatea g(x(T)) dx — dr = dr din care printr-o schimbare de variabila se obtine Г / —— du = t + C Jx* g(u) ( ) Rezulta in acest fel ca o soluție x = x(t) a ecuației diferentiale ( ) este solutie pentru ecuatia implicita G(t, x; C) = ( ) in care x G(t,x; C) = t + C — / du ( - ) Jx* g(u) Este usor de aratat folosind teorema functiilor implicite ca, daca x(t; C) este o solutie a ecuatiei ( ), atunci este si solutie a ecuatiei diferentiale ( ) Observația Daca functia g se anuleaza în x* G J, atunci functia constanta x(t) = x* este solutie a ecuatiei diferentiale ( ) Observația Pentru to G IR si xo G J, problema determinarii acelei solutii a ecuatiei ( ) care verifica conditia suplimentara x(to) = xo se numeste problema Cauchy sau problema cu date initiale: x = g(x) x(to) = xo ( ) iar solutia acesteia, x = x(t; to,xo), este data de ecuatia implicită: du — t — to ( ) Intr-o problema Cauchy to si xo sunt considerate cunoscute si se numesc date initiale Ecuații diferențiale autonome x = g(x) Concluzii Exista probleme de fizica care conduc la ecuații diferențiale de forma x = g(x) în care g este o funcție reala continua definita pe un interval deschis (c, d) C IR si nu se anuleaza Oricare ar fi solutia x(t) a ecuatiei diferentiale x = g(x) si oricare ar fi x* G (c, d) exista o constanta scalara C astfel încât x(t) este solutia C ecuatiei implicite / —— du — t — C = si reciproc, o solutie a acestei A* g(u) ecuatii implicite este solutie pentru ecuatia diferentiala Oricare ar fi t G IR si x G (c, d) exista o functie unica x = x(t) definita pe un interval deschis I (care contine pe t ) si cu valori in (c, d) care este solutia problemei cu date initiale x = g(x), x(t ) = x Daca functia g se anuleaza intr-un punct x* G (c, d) atunci functia constanta x(t) = x* este solutie a ecuatiei diferentiale Exerciții Sa se determine solutiile urmatoarelor ecuatii diferentiale (cu calculatorul): a) x = + x , x G IR R: b) x = e-x, x G IR R: c) x = k • x, x > R: d) x = k • x, x R: x(t) = tan(t + C), t + C = ( k + ) • П x(t) = ln(t + C), t + C > x(t) = C • ekt, C > , t G IR x(t) = C • ekt, C — , x > R: +t = C • et — x c) x +t o ^—, t > , x G IR x — R: x — arctan x = ln t+ + C x x Ecuații diferențiale cu variabile separate Sa se rezolve următoarele probleme Cauchy și sa se reprezinte soluțiile (cu calculatorul): a) x = y/tx, t > , x > , t = , x = R: x(t) = / — X b) x = — t , tGIR , xG ( , ), to = , x = R: — t — ^ = c) — x x = , t — t x c) t * = - — x ' - / \ (t+ ) R: x(t)= ln(t+ ) + —+ C) • t+ —t d) x = x — t R: x(t) = t + t + et C CAPITOLUL Sa se rezolve următoarele probleme cu date inițiale și sa se reprezinte grafic solutiile lor (cu calculatorul): e , q ^ ~i a) x = — tx + t , t = , x = R: x(t) = -t -\ e— + ’ ’ v b) X = - x — ln t, t to = , Xo = R: x(t)= f— ln + ) t c) X = —x + e*, t = , x = R: x(t) = e* + e * d) x = —ax+bept, t = , x = e R: x(t) = (be(p+ ^t—b +p+A— v a+p —a t Ecuatia diferențială a lui Bernoulli Ecuatia diferențială a lui Bernoulli Ecuația diferențială a lui Bernoulli are forma x = A(t) x + B(t) xa ( - ) în care funcțiile A și B sunt funcții reale continue A, B : (a, b) IR si se considera cunoscute, a este un numar real diferit de si cunoscut, iar funcția necunoscuta x(t) este pozitiva Pentru a determina soluțiile x (pozitive) ale ecuației ( - ) se introduce o noua funcție necunoscuta y = x -a- Aceasta verifica ecuația: ■^ = ( — a) A(t) y + ( — a) B(t)- ( - ) Ecuația ( - ) este o ecuație liniara de ordinul întâi si soluțiile ei sunt date de formula: y(t) = C e( -a^‘*a(t )dT + ( - ) / c t \ + ( ( — a) B(u) e-( -a) t* A(T)dTdu) e( -a)^- A(T)dT \ Jt* / Soluțiile pozitive x(t) ale ecuației ( - ) se determina din y(t) cu formula x(t) = y(t) si in general sunt definite pe (a,b)- Pentru t G (a,b) si xo > ecuația ( - ) are o soluție care verifica x(to) = xo si este data de formula x(t; to,xo)= y (t; t ,x ) ( - ) unde: c t y(t; to,yo) = yo e( -aH* A(T)dT + ( — a) B(u) e-( -a)^A(T)dTdu ( - ) Jt si yo = xo -a - Observația Ecuația Bernoulli apare în studiul miscarii corpurilor în medii care opun o rezistența la miscare de forma R = k v + k va, v fiind viteza corpului- CAPITOLUL Problema Sa se determine curba r = r(u) știind ca aria sectoarelor limitate de curba, raza vectoare a punctului P (r , u ) si raza vectoare a punctului P(r, u) este proportionala cu produsul ru, coeficientul de proporționalitate fiind a Rezolvare: Conform enuntului avem: r du = aru din care prin derivare obtinem: r = a (r u + r) care este o ecuatie Bernoulli Concluzii Exista probleme de fizica care conduc la ecuatii diferentiale de forma x = A(t) x + B(t) xa, (a G IR , a = , ) (numita ecuatia diferentiala a lui Bernoulli) in care A, B sunt functii reale continue definite pe un interval I C IR O functie pozitiva x = x(t) este solutie a ecuatiei Bernoulli daca si numai daca functia y(t) = [x(t)] —* este solutie a ecuatiei diferentiale liniare de ordinul întâi y = ( -a) A(t) y + ( — a) B(t) Determinarea solutiilor pozitive ale ecuatiei Bernoulli se reduce la rezolvarea unei ecuatii diferentiale liniare de ordinul întai Ecuatia diferențială a lui Bernoulli Exerciții Sa se determine soluțiile pozitive ale ecuațiilor: a) у у | y si o functie (t , x ) si de rază r astfel încat dP dQ z și , atunci pentru orice (t ,x ) reală U = U(t,x) definită pe discul centrat să aibe loc relația ( ) in G Demonstrație: Pentru un punct (t ,x ) discul centrat în (t , x ) si de raza r > sa fie inclus în Q Pornind de la faptul dU />t ca pe disc trebuie sa avem —— = P deducem ca U(t,x)= P(t, x)dr+Ф(я) dt Jto dU unde Ф este o functie de clasa C necunoscuta Impunand conditia — = Q dx obtinem egalitatea: Q se considera r > astfel ca C dP / -г—(t,x) dr + Ф'Ы = Q(t,x) Jt dx Tinand seama acum de egalitatea deducem ca: dx dt [ -^-(r,x) dr + Ф'(ж) = Q(t,x) to dt Efectuand integrarea se obtine egalitatea: Q(t, x) - Q(t , x) + Ф'(ж) = Q(t, x) din care rezulta: Ф'(я) = Q(t ,x) Prin urmare functia Ф(я) este data de formula: Ф(я) = /* Q(t ,y) dy + C Jxo ( ) CAPITOLUL în care C este o constanta reala Revenind la funcția U(t,x) obținem ca aceasta este data de formula: f’t px U(t,x) = P(t, x) dr + / Q(t ,y) dy + C J to J Xq ( ) Formula aceasta defineste o multime de functii U(t, x) care au proprietatea exprimata prin relatia ( ) dP dQ Comentariu: Propozitia arata ca egalitatea - — = - — este o conditie dx dt suficienta pentru ca sa existe în vecinatatea oricarui punct (t ,x ) G Q o functie U(t, x) de clasa C astfel ca dU = P dt + Q dx Mentionam ca si reciproca acestei afirmatii este adevărata Mai precis este adevarata urmatoarea afirmare: daca exista r > si o functie U(t, x) de clasa C pe discul centrat în (t , x ) si raza r astfel ca dU = P dt + Q dx pentru orice (t,x) din acest disc, atunci functiile P si Q sunt de clasa C dP dQ si — = — pentru orice (t, x) din disc Acest rezultat se obtine folosind dx dt posibilitatea inversarii ordinii de derivare, stabilit de Schwartz Observația Daca functiile P si Q sunt de clasa C pe Q C IR dP dQ dar —— = atunci ecuatia ( ) nu este o ecuatie cu diferentiala totala dx dt exactăa si metoda prezentatăa nu poate fi utilizatăa pentru determinarea solutiilor ecuatiei In acest caz este util sa observam ca ecuatia ( ) are aceleasi solutii ca si ecuatia P(t, x) • y(t, x) x — Q(t, x) • y(t, x) în care y(t,x) este o functie de clasa C care nu se anuleaza ( ) Datorita acestui fapt apare natural sa încercam sa determinam functia y(t, x) astfel ca ecuatia ( ) să fie cu diferentiala totala Impunand aceasta conditie rezulta ca functia y(t,x) trebuie sa verifice relatia: dP , dQ + Q dx + P dx dt + Q dt ( ) O functie care verifica ( ) se numeste factor integrant, iar relatia de dependenta functionala ( ) se numeste ecuatia factorului integrant Ecuatii cu diferențială totală exactă Factor integrant Exerciții Sa se rezolve următoarele ecuații cu diferențiale totale: tx a) x = — R: t x(t) + etx(t) = C ! tetx - t v ! b) tm + tx + x = xn + t x+ x R: x(t)n+ m+ n+ +t x(t) +ln(tx(t)) = C tx — x „ , o „ c) x = — R: t x(t) — tx(t) = C t — tx Sa se rezolve urmatoarele probleme Cauchy: a) x = -І, to = , xo = R: x(t) = t + -/ t + t—x t b) x = -, t = , x = R: x(t) = - —t + x Sa se rezolve ecuatiile diferentiale stiind ca ele admit factor integrant h = ^(t): t sin x + x cos x t a) x = R: ^(t) = et t cos x — x sin x et [(t — ) sinx(t) + x(t) cos x(t)] = C b) — t x t (x — t) R: h(t) = x(t) — t x(t) — = C x = CAPITOLUL Sa se rezolve ecuațiile diferențiale știind ca ele admit factor integrant p = p(x): a) x = x( — t x) —t R: p(x)= x b) x = t x x — t + t x(t) t =C R: p(x) x t x(t) + x(t) = C Calculul simbolic al soluțiilor ecuațiilor de ordinul întâi Calculul simbolic al soluțiilor ecuațiilor diferențiale de ordinul întâi Foarte multe dintre modelele matematice ale unor fenomene din realitate conțin cel putin o ecuație diferențiala Toate softurile comerciale de matematica (Maple, Mathematica, Mathcad) ofera posibiltatea sa rezolvam numeric aceste probleme Exemplele de rezolvare numerica care sunt în acest curs vor fi prezentate în programul Maple , versiune care acopera toate celelate versiuni de Maple în momentul de fata Pentru rezolvarea numerica a ecuatiilor diferentiale cu programul Maple se foloseste functia dsolve (solve ordinary differential equations - ODEs) cu una din urmatoarele sintaxe : dsolve(ODE); dsolve(ODE, x(t), extra args); dsolve({ODE, ICs}, x(t), extra args); în care: ODE - ecuatia diferentiala ordinara pe care dorim sa o rezolvam x(t) - functia necunoscuta pe care dorim sa o determinam ICs - conditiile initiale extra args - argumente optionale care se folosesc pentru schimbarea formei de afisare a solutiei (explicita, implicita, parametrica), a metodei de rezolvare a ecuatiei (separarea variabilelor, Bernoulli, Riccati, etc ) Pentru exemplificare, consideram ecuatia diferentială de ordinul întai: x —— • ( - x); t G R -{- }, x G R -{ } ( ) +t Aceasta ecuatie este cu variabile separate (caz particular de ecuatie liniara) Prin utilizarea sintaxei dsolve(ODE) se obtine multimea solutiilor ecuatiei date (ecuatia familiei de curbe integrale scrisa sub forma explicita): CAPITOLUL > dsolve(diff(x(t),t)=(t/( +t))*( -x(t))); x (t) = + -C^) (e-t + e-tt) Daca dorim ca soluțiile sa fie afișate sub forma parametrica, atunci se folosește argumentul optional ‘parametrici si obtinem: > dsolve(diff(x(t),t)=(t/( +t))*( -x(t)),x(t),'parametric'); x (t) = - ~CÎ e-tt ai • Se mai poate utiliza ca argument optional ”metoda de rezolvare a ecuatiei” Daca dorim sa se rezolve ecuatia diferentiala ca o ecuatie liniara, atunci se foloseste argumentul optional [linear] si obtinem: > dsolve(diff(x(t),t)=(t/( +t))*( -x(t)),x(t),[linear]); x (t) = ( + + (e-t + e-tt), iar daca dorim sa se rezolve ecuatia diferentiala ca fiind o ecuatie cu variabile separate, atunci folosim argumentul optional [separable] si obtinem: > dsolve(diff(x(t),t)=(t/( +t))*( -x(t)),x(t),[separable]); x (t) = ( ai et- -t)e-t ai Nespecificand metoda de rezolvare Maple va alege una dintre ele Deoarece in secventele de mai sus nu s-a dat nici o conditie inițială, Maple a afisat raspunsul cu ajutorul unei constante necunoscute Daca specificam si conditia initiala atunci calculatorul va rezolva o problema cu conditii initiale (Problema Cauchy) si va afisa solutia acesteia Pentru ecuatia diferentiala ( ) vom considera doua Probleme Cauchy deoarece domeniul de definitie al membrului drept este reuniunea (- x , - ) x IR U (- , +x) x IR Daca consideram t > - si conditia initiala x( ) = , atunci se obtine solutia: > dsolve({diff(x(t),t)=(t/( +t))*( -x(t)),x( )= },x(t)); x (t) = ( el / e-e -^ (e-t + e-tt) x (t) — I +t / e~ ) (e + e ^), iar daca consideram t dsolve({diff(x(t),t)=(t/( +t))*( -x(t)),x(- )= },x(t)); Calculul simbolic al soluțiilor ecuațiilor de ordinul întâi x (t) = (i+ + e~ } (e-t + e-tt) Pentru reprezentarea grafica a soluției unei probleme cu date inițiale, programul Maple folosește funcția plot (create a two-dimensional plot of func-tions) Utilizarea acesteia implica următoarea sintaxa: plot(f,h,v); în care: f - functia care trebuie reprezentata grafic; h - domeniul de definire al functiei pe axa orizontala; v - (optional) domeniul de variatie al functiei pe axa verticala Solutia Problemei Cauchy a ecuatiei ( ) corespunzatoare conditiei initiale x( ) = este reprezentata pe Figura > f :=(exp(t)/( +t)- / *(exp(- )*exp( )- )/exp(- ))* (exp(-t)+exp(-t)*t): > plot(f ,t=- infinity); Figura iar solutia Problemei Cauchy corespunzatoare conditiei initiale x(- ) = este reprezentatăa pe Figura CAPITOLUL Figura Dupa cum se poate observa din instrucțiunile de mai sus s-a atribuit variabilei f functia solutie a Problemei Cauchy si apoi am folosit in instructiunea plot In general, este recomandabil sa se atribuie unor expresii matematice variabile, deoarece aceasta simplifica scrierea In cele ce urmeaza, continuam exemplificarea rezolvand trei probleme cu date initiale si, în fiecare caz, vom reprezenta grafic solutia: Ecuatia liniara X = — x + e* ( ) > dsolve(diff(x(t),t)=-x(t)+ *exp(t),x(t),[linear]); x (t) = e* + e-t-C > dsolve({diff(x(t),t)=-x(t)+ *exp(t),x( )= },x(t), [linear]); x (t) = e* + e-t > plot(exp(t)+exp(-t),t=- ); Calculul simbolic al soluțiilor ecuațiilor de ordinul întâi - - Figura sau > plot(exp(t)+exp(-t),t=- ,color=black,style=point, axes=boxed); Figura în care am folosit diferite comenzi opționale referitoare la modalitatea de afișare a graficului Ecuatia de tip Riccati x = —x + - • x t? ’ ( ) t > CAPITOLUL > eq:=diff(x(t),t)=-x(t)" +( /t)*x(t)- /t" ; — L,- (t) — C* (t)) + x(t) t- eq ♦ Ltx (t) (x (t)) I t t > dsolve(eq,'explicit',[Riccati]); x (t) = ( C - / t- )- / + / > dsolve(eq,[Riccati]); x (t) = ( C - / t- )- t- + t- > dsolve({eq,x( )= },x(t)); x (/) = t + X (/) ( +t )t > dsolve({eq,x( )= },x(t),[Riccati]); x (t) = (- / - / - )- - + - > sol :=( *t" + )/(( +t" )*t): > sol := /((- / - / /t~ )*t" )+ /t: > plot([sol ,sol ],t= ,x= ,color=[red,blue], style=[point,line]); Figura Calculul simbolic al soluțiilor ecuațiilor de ordinul întâi în secvențele de mai sus observam ca, argumentul opțional in care cerem sa se afiseze soluția sub forma explicita este inutil, deoarece acest lucru este făcut automat de dsolve Deasemenea, daca folosim argumentul optional [Riccati] solutia ecuatiei difera doar aparent (cele doua solutii afisate coincid dupa cum se poate observa din Figura unde am reprezentat simultan ”ambele” forme ale solutiei in acelasi sistem de coordonate) Ecuatia cu factor integrant x = - x V+ ’ x - t + = ( - ) > dsolve(diff(x(t),t)=- *t*x(t)/( *x(t)~ -t~ + ),‘explicit‘; x (t) = - / C ± / у/ C - t + > dsolve(diff(x(t),t)=- *t*x(t)/( *x(t)~ -t~ + ),‘implicit‘); + x (t) — (x (t)) + -C = > dsolve({diff(x(t),t)=- *t*x(t)/( *x(t)~ -t~ + ),x( )= }) ; x (t) = / V — t > plot( / *( - *t~ )~( / ), t=- ); Figura CAPITOLUL Ecuația cu factor integrant a fost rezolvata de Maple fara specificarea factorului integrant p = p(t,x) iar solutia a fost afisata sub forma explicita in primul caz, respectiv sub formă implicită în al doilea caz In Figura este reprezentată solutia Problemei Cauchy corespunzatoare Capitolul Ecuații diferențiale de ordin superior rezolvabile prin metode elementare Definiția O ecuație diferențială de ordinul n > este o relație de dependență funcțională de forma g(t,x,x, , x(ra)) = ( ) înțre funcția idențica t t definița pe un ințerval I C IR necunoscuț, o funcție necunoscuța x(t) si derivațele ei x ,x(ra) pană la ordinul n definițe pe acelasi ințerval în ecuația ( ) funcția g se considera cunoscuta și rezolvarea ecuației înseamnă determinarea funcțiilor necunoscute x care verifica ecuația Definiția O funcție reală x de clasă Cn definiță pe un ințerval deschis I C IR se numește soluție a ecuației ( ) dacă pentru orice t G I, sisțemul ordonaț (t,x(t),x(t), ,x(n\t)) aparține domeniului de definiție a lui g si g(t, x(t),x(t), , x(ra)(t)) = ( ) Vom prezența cațeva cazuri de asemenea ecuații care se rezolva cu mețode elemențare si probleme concrețe din diferițe domenii care au condus la asemenea ecuații CAPITOLUL Ecuații diferențiale liniare de ordinul al doilea cu coeficienți constanți Problema Sa se determine variația curentului intr-un circuit format dintr-o rezintenta R, o bobina cu inductanta L si un condensator de capacitate C legati în serie si conectafi la o sursa de curent alternativ de tensiune electromotoare E = E • cos xl Rezolvare: Fie i(t) intensitatea curentului din circuit la momentul t Caderile de tensiune pe elementele circuitului sunt: uR = R • i; di UL =L • dt; UC = C J i(t)dt; conform celei de-a doua legi a lui Kirchhoff, suma caderilor de tensiune pe bobina, rezistenta si condensator este egala in orice moment cu tensiunea electromotoare a generatorului Prin urmare avem: di f L • — + R • i + — / i(t)dt = E • cosші, dt j iar prin derivare se obtine ca intensitatea a curentului verifica egalitatea: T d i di L • dt +R • dt+ C = -Eo •ш •'■ w Prin urmare avem de determinat o functie i(t) care împreuna cu derivatele ei de ordinul întâi si doi verifica relatia de dependenta functionala (*) In (*) cu exceptia functiei i totul este cunoscut Pentru a determina functia necunoscuta i(t) vom arata în continuare cum se rezolva o ecuatie diferentiala liniara de ordinul al doilea cu coeficienti constanti O ecuatie diferentiala liniară de ordinul al doilea cu coeficienti constanti este o ecuatie diferentiala de forma: a X + a X + a x = f (t) ( ) în care a ,a , a sunt constante reale cunoscute, a = , f (t) functie continua cunoscuta si x este o functie reala de clasa C necunoscuta Ecuatii diferențiale de ordinul al doilea cu coeficienți constanți Observația Daca f = atunci ecuația ( ) se numește ecuație diferențială liniara de ordinul doi cu coeficienți constanți omogena, iar daca f = ecuația ( ) se numesțe ecuație diferențiala liniară de ordinul doi cu coeficienți consțanți neomogena Vom dețermina mai înțai soluțiile ecuației omogene urmand apoi sa de-țerminam si soluțiile ecuației neomogene Fie ecuația omogena ațasața ecuației ( ): a x + a x + a x = ( ) Daca a = ațunci ecuația ( ) esțe o ecuație liniara de ordinul înțai: a X + a x = si soluțiile ei sunț dațe de formula /X a t x(t) = Ce-t în care C esțe o consțanța reală oarecare Observam ca raporțul — din ai exponenț, esțe soluția ecuației algebrice a • A + a = , iar la formula soluției x(t) = Ce ai se poațe ajunge nu numai pe calea descrisa în Capitolul § ci si cauțand soluții de forma x(t) = CeXt Aceasța esțe ideea pe care o vom folosi pențru a dețermina soluțiile ecuației ( ) Impunand unei funcții de forma x(t) = CeXt sa verifice ecuația ( ) rezulța ca A țrebuie sa verifice ecuația de gradul al doilea: a A + aiA + aox = ( ) Daca radacinile A si A ale ecuației ( ) sunț reale si disțincțe, ațunci funcțiile x (t) = C eX t si x (t) = C eX t sunț soluții ale ecuației ( ) si funcția x(t) = CieX t + C eĂ t esțe de asemenea soluție a ecuației ( ) Mai mulț, pențru orice t , x , x G R puțem dețermina în mod unic consțanțele C si C asțfel încaț sa aiba loc x(t ) = x si x(t ) = x ( ) CAPITOLUL în adevăr, impunând condițiile ( ) funcției x(t) = C eX t + C eX t, obținem următorul sistem de ecuații algebrice: = CieX t + C eX t = CieĂ t + C eX t° în care necunoscutele sunț Ci si C Determinantul acesțui sisțem esțe e(X +X )t • (A — Ai) si esțe nenul (Ai = Ă ), fapt pentru care sistemul are o soluție unica în particular rezulta de aici ca formula: x(t) = CieX t + C eX t ( ) reprezinta toate soluțiile ecuației ( ) în cazul în care ecuația ( ) are raadaacini reale distincte Daca ecuația ( ) are radacinile confundate Ai = A = A atunci pe langa funcția xi(t) = CieXt si funcția x (t) = C t • eXt este soluție a ecuației ( ) Prin urmare orice funcția x(t) de forma x(t) = CieXt + C t • eXt adicaa x(t) = eXt • (Ci + C t) ( ) este soluție a ecuației ( ) Mai mult, pentru orice t ,x ,x G IR putem determina în mod unic constantele Ci si C astfel încat sa aibe loc x(t ) = x si x(t ) = xj în adevar, impunand aceste condiții funcției data de ( ) rezulta urmatorul sistem de ecuatii algebrice: x = eXt (Ci + C to) x = AeXt • Ci + C eXt + eXt • to • A • C al carui determinant este e Xt = în particular rezultaa de aici caa, formula ( ) reprezintaa toate solutiile ecuației ( ) în cazul în care ecuația ( ) are radacinile confundate Ramane sa consideram cazul în care ecuația ( ) are radacinile complex conjugate Ai = + iv si A = — iv în acest caz consideram funcțiile xi(t) = CieMt • cos vt si x (t) = C eMt • sin vt Ecuatii diferențiale de ordinul al doilea cu coeficienți constanți (Ci, C constante reale) și arătăm că fiecare din acestea este soluție a ecuației ( ) Demonstratia se face prin verificare Pentru exemplificare facem acest calcul in cazul functiei xi(t): x (t) = Ci^ • • cos vt — Civ • • sin vt xi(t) = Ci^ • • cos vt — Ci^v • • sin vt — Civ • • cos vt si înlocuind în ecuatia ( ) avem: a x + aiâ: + a xi = Ci • • cos vt [(p — v )a + ^ai + a ^ + + C • • sin vt [— p,va — vai] Deoarece a (^ + iv ) + a (p, + iv) + a = avem (^ — v )a + ^a- + a + i [ ^va + vai] = si prin urmare: (^ — v )a + ^a- + a = si ^va + vai = Ținand seama de aceste egalitati deducem egalitatea a Xi + aix i + aoxi = care arata ca functia xi (t) = Ci • • cos vt este solutie a ecuatiei diferentiale ( ) La fel se arata că functia x (t) = C • • sin vt este solutie a ecuatiei diferentiale ( ) Astfel, rezulta ca orice functie x(t) = Ci • • cos vt + C • • sin vt ( ) este solutie a ecuatiei ( ) Aratam în continuare ca pentru orice to,x ,x G IRI putem determina constantele Ci si C în mod unic astfel încat sa aibe loc x(t ) = x si x(t ) = xi CAPITOLUL Impunând aceste condiții funcției ( ) rezulta următorul sistem de ecuații algebrice: = eMt° • [Ci • cos vt + C • sin vt ] = eMt° • [C • (^ cos vt — v sin vt ) + C • (^ sin vt + v cos vt )] avand ca necunoscute constantele C , C Determinantul acestui sistem algebric este v • e t° si este diferit de zero In particular, rezulta de aici ca formula ( ) reprezinta toate solutiile ecuatiei ( ) in cazul în care ecuatia ( ) are radacinile complexe Am ajuns in acest fel sa determinam toate solutiile ecuatiei ( ) Aceasta însa nu permite înca sa rezolvam problema pusa la începutul paragrafului, pentru ca aceasta conduce de fapt la ecuatia ( ), adica: a x + a x + a x = f (t) în care functia f este data Reamintim ca, deosebirea dintre ecuatiile ( ) si ( ) consta în faptul ca în membrul drept al ecuatiei ( ) este o functie continua care nu neaparat este functia identic nula, adica este o ecuatie diferentiala de ordinul al doilea cu coeficienti constanti neomogena Pentru determinarea solutiilor ecuatiei ( ) este important sa observam la început ca, daca x(t) este o solutie fixata a ecuatiei ( ) si x(t) este o solutie oarecare a aceleiasi ecuatii, atunci diferenta x(t) = x(t) — x(t) este o solutie oarecare a ecuatiei ( ) Intrucat solutiile x(t) ale ecuatiei ( ) sunt cunoscute, determinarea solutiilor x(t) ale ecuatiei ( ) revine la determinarea unei singure solutii x(t) ale acestei ecuatii O solutie particulara x(t) pentru ecuatia ( ) se determina cu metoda variatiei constantelor a lui Lagrange (un procedeu asemanator cu cel descris în Cap § ) In continuare prezentam aceasta metoda în cazul în care ecuatia algebrica ( ) are radacinile reale distincte Л , Ă In acest caz solutiile ecuatiei omogene ( ) se scriu sub forma ( ): T(t) = CieX t + C eĂ t Ecuatii diferențiale de ordinul al doilea cu coeficienți constanți Soluția particulara x(t) a ecuației neomogene ( ) se cauta sub aceeași forma considerând însă C , C funcții de clasa C de variabila t: x(t) = Ci(t)eX t + C (t)eX t ( ) Pentru a impune funcției x(t) sa verifice ecuația ( ) calculam derivata acesteia si obținem: X(t) = Ci(t)eX t + C (t)eX t + Ci(t)XieX t + C (t)X eX t ( ) în continuare ar trebui sa calculăm derivata de ordinul al doilea x prin derivare în raport cu t în expresia ( ) Aceasta ar introduce derivatele de ordinul al doilea ale funcțiilor C (t),C (t) de existența carora nu ne-am asigurat De aceea impunem condiția suplimentara: Ci(t)eX t + (C (t)eX t = ( ) Cu aceasta ( ) devine: x(t) = Ci(t)XieX t + C (t)Ă eX t ( ) iar prin derivare obținem: x(t) = Ci(t)XieX t + C (t)X eX t + Ci(t)X?eX t + C (t)X eX t ( ) înlocuind ( ) si ( ) în ( ) rezulta: Ci (t)(« X + aiXi + ao)eX t + C* (t)(a X + aiX + ao)eX t+ + C i(t)a XieX t + (C (t)a X eX t = f (t) C i(t)XieX t + (C (t)X eX t = f (t) ( ) a Astfel, sistemul de ecuații algebrice format din ecuațiile ( ) si ( ): Ci(t)eX t + C (t)eX t = i ( ) CJi(t)XieX t + CJ (t)X eX t = - f(t) CAPITOLUL în care necunoscutele sunt C (t) (derivatele funcțiilor C (t) și C (t)), are determinantul (A — A )e(X +X )t = si permite determinarea funcțiilor CCi(t) si C (t): C (t) = = (A A ) • e-(X +X )t • eX t • f(t) a (A — A ) ( ) CC (t) = = A ) • e- ) = r x(n )(tn) = xn x( o ) = xo , x( o) = xo , , x ( o) = xo In particular rezulta de aici ca formula ( ) reprezinta toate soluțiile ecuației ( ) în acest caz CAPITOLUL Daca printre rădăcinile ecuației caracteristice ( ) există și rădăcini complexe simple, de exemplu A = p + iv si A = p — iv, atunci fiecarei perechi de radacini complex conjugate îi corespund solutiile x\ (t) = Cx • • cos vt si x x (t) = Cx • eMt • sin vt Pentru p = aceste solutii devin: xX(t) = Cx • cos vt si xX(t) = Cx • sin vt Astfel, daca ecuatia caracteristica are k radacini complexe simple Aj = Pj + ivj si Aj = pj — ivj, j = , k si n — k radacini reale simple A k+ , , An, atunci orice functie x(t) data de: x(t) = ^ Cj • e^jj • cos vjt + ^ Cj • e^jj • sin vjt + Cj • eXj ( ) j= j= j= k+ este solutie a ecuatiei ( ) (Cj,Cj ,j = ,k si Cj, j = k + ,n sunt constante reale oarecare) Mai mult, oricare ar fi t ,x ,x , ,хП- G IR putem determina in mod unic constantele Cj , Cj, j = , k si Cj, j = k + , n astfel încât sa aiba loc x(t ) = x , x(t ) = x , , x(n- )(t ) = x£- în particular rezulta de aici ca formula ( ) reprezinta toate solutiile ecuatiei ( ) în acest caz Daca ecuatia caracteristica ( ) are k radacini reale A , , Ak avand ordine de multiplicitate q , , qk si l radacini complex conjugate p ±iv , , pl± ivi avand ordine de multiplicitate r , ,ri, atunci orice functie x(t) data de formula: k i x(t) У eXjj • Pqj- (t) + У eMj* • [Qr,- (t) • cos vjt + Rrj- (t) • sin vjt] j= j= ( ) este solutie a ecuatiei ( ), unde Pqj- (t) sunt polinoame de grad qj — cu coeficineti reali nedeterminati si Qrj- , Rrj- sunt polinoame de grad rj — cu coeficienti reali nedeterminati Mai mult, oricare ar fi t ,x ,xj, ,x((- G IR putem determina în mod unic coeficientii polinoamelor Pqj- , Qqj- , Rqj- astfel încat sa aiba loc x(t ) = x , x(t ) = xj, , x(n- )(t ) = x£- în particular rezulta de aici ca formula ( ) reprezinta toate solutiile ecuatiei ( ) în acest caz Ecuatii diferențiale liniare de ordinul n cu coeficienți constanți Reamintim ca obiectul acestui paragraf este rezolvarea ecuației diferențiale de ordinul n, cu coeficienti constanti neomogena ( ): arax(ra) + ara x(ra ) + + a^r + aox = f (t) în care a , a , , an , an sunt constante reale date, an = , f functie cunoscuta continua si x este functie reala de clasa Cn necunoscuta Pentru determinarea solutiilor ecuatiei ( ) este important sa observam ca daca x(t) este o solutie fixata a ecuatiei ( ) si x(t) este o solutie oarecare a aceleiasi ecuatii, atunci diferenta X(t) — x(t) = x(t) este o solutie oarecare a ecuatiei diferentiale liniare omogene cu coeficienti constanti, ( ) Intrucat solutiile x(t) ale ecuatiei omogene ( ) sunt date in general de ( ), determinarea solutiilor x(t) ale ecuatiei ( ) revine la determinarea unei singure solutii x(t) ale acestei ecuatii O solutie particulara x(t) pentru ecuatia ( ) se determina cu metoda variatiei constantelor a lui Lagrange, un procedeu asemenator cu cel descris în paragraful precedent Vom ilustra acest procedeu pe un exemplu (n = ): Exemplul Sa se determine solutiile ecuatiei: x + x + x = et Consideram ecuatia omogena x + x + x = Ecuatia caracteristica asociata este: A + A + A = ale carei radacini sunt: Ao = , Ai = — — i, A = — + i Solutiile ecuatiei omogene sunt date de: y(t) = C + C • e t • cos t + C • e t • sin t Cautam x(t), o solutie particulara pentru ecuatia neomogena, sub forma x(t) = C (t) + C (t) • e t • cos t + C (t) • e t • sin t CAPITOLUL în care C'i(t), C^t), Сз(і) sunt funcții de clasă C care trebuiesc determinate Calculăm derivata întâi a funcției x(t) și obținem: x = C*i + C • e~ t • cos t + C • e~ t • sin t — C ■ e~ t ■ cos t— — C ■ e~ t • sin t — С* • e~ t • sin t + C • e~ t ■ cos t Impunem ca Ci, C , C să verifice: C*i + C • e~ t ■ cos t + C • e~ t • sin t = și obținem: x = —C • e~ t ■ ( cos t + sin t) + C ■ e~ t ■ (cos t — sin t) Calculăm derivata a doua a funcției x(t) și obținem: x = — C • e~ t • ( cost + sini) + C • e~ t • (cost — sint)+ + C ■ e~ t ■ ( cos t + sin t) — C • e~ t • (cos t — sin t) — —C ■ e~ t ■ (— sint + cost) + C ■ e~ t • (— sint — cost) = = — C ■ e~ t ■ ( cost + sint) + C ■ e~ t ■ (cost — sint) + +C ■ e~ t ■ ( cost + sint) + C ■ e~ t ■ ( sint — cost) Impunem ca C , C să verifice: — C ■ e~ t ■ ( cost + sint) + C ■ e~ t ■ (cost — sint) = și obținem x = C • e~ t ■ ( cos t + sint) + C • e~ t ■ ( sint — cos t) De aici calculăm derivata a treia a funcției x(t) și obținem: x = C • e~ t • ( cost + sint) + C • e~ t • ( sint — cost) — — C ■ e~ t ■ ( cost + sint) — C ■ e~ t ■ ( sint — cost) + +C* • e~ t ■ (— sin t + cos t) + C • e~ t ■ ( cos t + sin t) = = C ■ e~ t ■ ( cost + sint) + C ■ e~ t ■ ( sint — cost) + +C ■ e~ t ■ (— cost — sint) + C ■ e~ t ■ (— sint + cost) Ecuatii diferențiale liniare de ordinul n cu coeficienți constanți înlocuind toate acestea în ecuația data rezulta: C • e- t • ( cost + sint) +C • e- t • (- cos t — sin t) +C • e- t • ( cos t + sin t) —C • e- t • ( cos t + sin t) +C • e- t • ( sint — cos t) +C • e- t • (— sin t + cos t) +C • e- t • ( sin t — cos t) +C • e- t • ( cos t — sin t) e* + + sau C • e t • ( cost + sint) + C • e t • ( sint — cost) = e* Aceasta egalitate împreuna cu sistemul de conditii impus pe parcurs functiilor C ,C ,C conduce la urmatorul sistem liniar de ecuatii algebrice în necunoscutele C , C , C : C +C • e t • cos t +C • e t • sin t = , avem ca x(t) = y(ln t) iar prin derivare succesiva obținem: dx dy dt dT t d x d y dy f d y dy\ dt dT t dT t t \dT dTj d x f d y dy\ (d y d y\ dt t \d,T dTjrt \dT dT ) = A d y d y + țy\ t \dT dT dT) Daca presupunem ca pentru : x(t) = y(ln t)- Pentru t ) Maple folosește aceeași funcție dsolve (solve ordinary differential equations - ODEs) care a fost prezentata în capitolul anterior Noutatea care apare aici consta în scrierea sintaxei pentru derivatele de ordin superior De exemplu, derivata de ordinul al doilea a functiei x(t) poate fi scrisa într-unul din urmatoarele moduri: diff(x(t),t,t) diff(x(t),t$ ) (D@@ ) (x) (t) Pentru exemplificare, vom rezolva cateva ecuatii si probleme cu date initiale: Ecuatia diferentiala liniara de ordinul al doilea cu coeficienti constanti omogena: X — X + x = ; ( ) > > > > > eq :=diff(x(t),t,t)- *diff(x(t),t)+ *x(t)= ; eq := x (t) — dx (t) + x (t) = dsolve(eq ,x(t)); x (t) = C e + e dsolve({eql,x( )= ,D(x)( )= },x(t)); x (t) = — e sol := *exp( *t)/exp( )- *exp( *t)*t/exp( ): plot(sol ,t=-infinity infinity); CAPITOLUL Se observa ca, daca nu s-a dat nici o condiție inițială soluția generala este afișata cu ajutorul a doua constante Mai precis, numarul constantelor este acelasi cu ordinul ecuatiei, in cazul ecuatiei de ordinul al doilea solutia generala exprimandu-se cu ajutorul a doua constante Pentru ca Maple sa afiseze solutia unei Probleme Cauchy în cazul unei ecuatii de ordin superior trebuie sa-i dam n conditii initiale: x(t>) — x xft>) — x ) — xn- Ecuatia diferentiala liniara de ordinul al patrulea cu coeficienti constanti omogena: x( ) - x + x — ; ( ) > eq :=diff(x(t)(x(t),t,t)+ *x(t)= ; eq :— ^ x (t) - x (t) + x (t) — > eq :=diff(x(t),t$ )- *diff(x(t),t$ )+ *x(t)= ; eq :— ^ x (t) - x (t) + x (t) — > eq :=(D@@ )(x)(t)- *(D@@ )(x)(t)+ *x(t)= ; eq :— (D( )) (x) (t) - (D( )) (x) (t) + x (t) — > dsolve(eq ,x(t)); x (t) — C - + C e-t + C e + C/ С > dsolve({eq ,x( )= ,D(x) ( )= ,(D@@ ) (x) ( )= , (D@@ )(x)( )= },x(t)); x (t) — - / e-t + / e- + / e - / e Calculul simbolic al soluțiilor ecuațiilor de ordinul n > sol :=- / *exp(-t)+ / *exp(- *t)+ / *exp( *t)- / *exp(t): > plot(sol ,t=- ); pentru o ecuație diferențială de ordinul patru s-au folosit patru condiții inițiale Se observa deasemenea, sintaxa corespunzatoare derivatelor de ordin superior a fost scrisa în cele trei moduri prezentate la începutul paragrafului Ecuatia diferentiala liniara de ordinul al treilea cu coeficienti constanti neomogena: x — x + x — x = sin t; ( ) > eq :=diff(x(t),t,t,t)- *diff(x(t),t,t)+ *diff(x(t),t) - *x(t)=sin(t); eq := fsx (t) — x (t) + dx (t) — x (t) = sin (t) > dsolve(eq ); x (t) — Л cos (t) sin (t) + + e tt x (t) — — cos (t) — І sin (t) + e + e t + ud e t > dsolve({eq ,x( )= ,D(x)( )= ,(D@@ )(x)( )= }); x (t) = — - cos (t) — sin (t) + - e * + e tt — e tt x (t) cos (t) sln (t) + e + e L e L > sol :=- / *cos(t)- / *sin(t)+ / *exp( *t)+ / *exp( *t)*t- / *exp( *t)*t~ : > plot(sol ,t=- ); CAPITOLUL Figura în cele ce urmeaza, vom mai prezenta o alta funcție de plotare DEtools [DEplot] (plot Solutions to an equation or a system of DEs) pentru a vizualiza solutia acestei probleme cu date initiale Utilizarea acesteia nu necesita rezolvarea ecuatiei in avans deoarece ecuatia este inclusa direct în instructiune Sintaxa acestei funcții poate avea una din urmatoarele forme: with(DEtools):DEplot(deqns, vars, trange, options); with(DEtools):DEplot(deqns, vars, trange, inits, options); with(DEtools):DEplot(deqns, vars, trange, xrange, yrange, options); with(DEtools):DEplot(deqns, vars, trange, inits, xrange, yrange, op-tions); in care: Calculul simbolic al soluțiilor ecuațiilor de ordinul n deqns - ecuația diferențială de orice ordin pe care dorim să o rezolvăm sau lista de ecuații diferențiale de ordinul întâi (in cazul sistemelor) vars - variabila independentă sau lista variabilelor independente trange - domeniul de definiție al variabilei independente inits - lista de condiții inițiale xrange - domeniul de variație al primei variabile dependente yrange - domeniul de variație al celei de-a doua variabile dependente options - diferite opțiuni: modul de afișare al soluției, metoda de rezolvare, etc > with(DEtools):DEplot(eq ,x(t),t=- ,[[x( )= ,D(x)( )= , (D@@ )(x)( )= ]],x=- ,stepsize= ,title='Solutia Problemei Cauchy'); Х(Т) ■ - - - t :- Figura Ecuația diferențială liniară de ordinul al treilea cu coeficienți variabili de tip Euler: t x + tx + x = Ini, i> ( ) > eq :=t~ *diff(x(t),t,t,t)+ *t*diff(x(t),t,t) + *diff(x(t),t)=ln(t); eq := і ^-ж (t) + t^x (i) + ftx (i) = ln (t) > dsolve({eq ,x( )= ,D(x)( )=l/ ,(D@@ )(x)( )= }); CAPITOLUL x(t) = - ln( ) + + (- + У ))-ln(t) + in( )- (in( )) - + i /^t(in(t)) - > sol :=- *ln( )+ +(- + *ln( ))*ln(t)/t+( *ln( )- *ln( )~ - )/t+ / *t*(ln(t)- ): > plot(sol ,t= infinity,axes=boxed); Figura Capitolul Sisteme de ecuatii diferențiale de ordinul întai liniare cu coeficienți constanti Sisteme de ecuatii diferentiale de ordinul întai liniare cu coeficienți constanți omogene Definitia Un sistem de ecuații diferențiale de ordinul întâi liniare cu coeficienți constanți omogen este un sistem de n relații de dependență funcțională de forma: Xi — an • xi + « • X + + ain • xn X — « • X + « • X + + « n ' Xn ( ) Xn — «ni • X + «n • X + + a nn • xn dintre un sistem de n funcții necunoscute x ,x , ,xn si derivatele acestora X ,X , ,xn In sistemul ( ) coeficienții sunt constante considerate cunoscute Definitia Un sistem ordonat de n funcții reale x ,x , ,xn de clasa CAPITOLUL C este soluție a sistemului ( ) dacă verifică dxi = ^ aij • xi(t) j=i pentru orice t G IR Definiția Fiind date t G IR' si (x°, x , , x^) G R problema determinării soluției (x (t), x (t), , xn(t)) a sistemului ( ) care verifica xi(t ') = x i = ,n se numeste problema cu date initiale sau problema Cauchy Pentru reprezentarea matriceala a sistemului ( ) notam cu A matricea patratica care are ca elemente constantele aij: A = (aijsi cu X matricea coloana X = (x , x , xn)T Cu aceste matrice sistemul ( ) se scrie sub forma matriceala: X = A • X ( ) în aceasta problema derivarea funcției matriceale X = X (t) înseamna derivarea elementelor matricei, iar produsul A • X însemna produsul dintre matricea A si matricea X Problema Cauchy (problema cu date initiale) se scrie matriceal sub forma: X = A • X, X(to) = X ( ) si consta în determinarea functiei matriceale X = X(t) care verifica ecuatia ( ) si conditia inițială X(t ) = X Teorema (de existentă a soluției problemei Cauchy) Pentru orice t G IR si X = (x ,x° , x°n)T problema Cauchy ( ) are o solutie definită pe IR Demonstrație: Consideram sirul de functii matriceale definite astfel: Sisteme de ecuații diferențiale de ordinul întâi liniare omogene X (t) = X = I • X X (t) = X + [ A • X (t)dr = I + (t - t ) A] • X X (t) = X + ! A • X (t)dT = [I + X (t) = X + j A • X (t)dT = = [/ I (t - t ) • A I (t - t ) • A = [ ' ! ' ! ! (t - t ) • A + (t - t ) • AX ! ! (t - t ) • A ] X ! ] + Xm(t) = X AXm- (r)dT = = [I I (t - t ) • A I (t - t ) • A = [ ' ! ' ! + + (t - t )m • Am] X m •X Funcțiile din acest șir verifica inegalitatea m+p lt — t |k II A IIk ||Xm+p(t) -Xm(t)|| A m= Sisteme de ecuații diferențiale de ordinul întâi liniare omogene Observația Funcția matriceala e(t-to)'A se numește matricea rezolvanta a sistemului ( ) O solutie a problemei Cauchy ( ) se obtine înmultind matricea e(t-to)'A cu matricea X : X(t; t ,X ) = e(t-to• X Aceasta solutie este definita pe IR Teorema (de unicitate a soluției problemei Cauchy) Problema Cauchy ( ) are o singura soluție Demonstrație: Presupunem prin absurd ca problema Cauchy ( ), pe langa solutia X(t; t , X ) determinata în teorema precedenta mai are o solutie X (t) Pe intervalul I de definitie a acestei solutii (I Э t ) scriem egalitatile : X(t; t ,X )= X + A ■[ X(t; t ,X )dr, (V)t G I dto si X(t) = X + A • /* X(t)dT, (V)t G I dto si deducem succesiv: X(t; t ,X ) - XC(t) = A • [\x(t; t ,X ) - XC(t)]dT Jto \\X(t; t ,X ) - XC(t)H T > (V)t G I Pentru t > t rezulta în continuare: ||X(t; t ,X ) - X(t)H - „ c -Z ° f + / I|A||-||X(t; t ,X ) - X'(T)BdT Jto CAPITOLUL ln ln l|A|| • ||X (t; to,X ) - X (t)H + f||A||J|X (т; to,X ) - X(r)||dr Jt U, > to =^ ||X(t; t ,X ) - X(t)H L, (V) > Pentru t fixat și rezulta ||X(t; t ,X ) - X(t)H = Astfel am aratat ca pentru orice t > t și t G I avem XC(t) = X(t; t , X ) Raționam analog pentru t si a + d , atunci orice soluție nenula a sistemului tinde la ( , ) ii) daca (a — d) + • b • c > si a + d > si a • d — b • c > , atunci orice soluție nenula a sistemului tinde in norma la +to iii) daca (a — d) + • b • c , atunci toate soluțiile nenule ale sistemului tind in norma la +to v) daca (a — d) + • b • c se obtine ca i(t) — i(t) pentru t ■ x adica, i(t) i(t) CAPITOLUL Exerciții Rezolvați următoarele ecuații diferențiale de ordin superior liniare cu coeficienți constanți prin metoda reducerii la un sistem de ecuații diferențiale de ordinul întai liniare cu coeficienți constanți: a) x — x = x( ) = x b) x + x + x = x( ) = x c) x — x + x = x( ) = x (П \ / = % e) x + x + x = x( ) = x ( ) = R: x(t) = e- + e - ( ) = R: x(t) = t • e ( ) = R: x(t) = e — te (П \ J= R: x(t) = sin t ( ) = R: x(t) = e- - sin( ^ ) a) x — x — x + x = x( ) ab( ) = x( ) = R: x(t) = • e- + • e - — • e - b) x — x + x — x = x( ) = x( ) = x( ) = R: x(t) = e- + (sin ) • sint — (cos ) • cos t c) x( ) — x + x = x( ) = ab( ) = x( ) = x ( ) = R: x(t) = — • e- + • e - — j • e - + j • e - Calculul simbolic al soluțiilor sistemelor de ecuatii liniare Calculul simbolic al soluțiilor sistemelor de ecuații diferențiale de ordinul întâi liniare cu coeficienți constanți Pentru rezolvarea numerica a sistemelor de ecuații diferențiale de ordin întâi Maple folosește funcția dsolve (solve ordinary differential equations -ODEs) care a fost prezenta în capitolele precedente In scrierea sintaxei ecuatia diferentiala va fi înlocuita cu lista de ecuatii diferentiale de ordinul întai care formeaza sistemul de ecuatii, respectiv conditia initiala va fi înlocuita cu lista conditiilor initiale xi(t ) = x corespunzatoare fiecarei functii necunoscute xi(t), i = ,n: dsolve({ODE , ODE , ,ODEn}); dsolve({ODE , ODE ,ODEn}, {x (t), x (t), ,xn(t)}, extra args); dsolve({ODE , ODE ,ODEn, x (t )=x ,x (to)=x , ,хп(і )=хП}, {x (t), x (t), , xn(t)}, extra args); Pentru exemplificare, vom rezolva urmatoarele siteme de ecuatii diferentiale de ordinul întai liniare cu coeficienti constanti: Sistemul de doua ecuatii diferentiale liniare cu coeficienti constanti omogen: x = —x + x x = xi+ x ( ) Pentru acest sistem vom considera conditiile initiale x ( ) = , x ( ) = si solutia problemei cu date initiale va fi reprezentată grafic utilizand trei instructiuni de plotare: > sys Eq :=diff(x (t),t)=-x (t)+ *x (t); sys Eq := jx (t) = —x (t) + x (t) > sys Eq :=diff(x (t),t)=x (t)+x (t); sys -Eq := jx (t) = x (t) + x (t) CAPITOLUL > dsolve(sysl Eql,sysl Eq ,xl( )= ,x ( )= ,xl(t),x (t)); {x (t) = / e t + / e t, xl (t) = / e t - / e~ t} > sol xl:= / *exp( *t)- / *exp(- *t): > sol x := / *exp( *t)+l/ *exp(- *t): > plot([sol xl,sol x ],t= ,color=[red,blue], style=[line,point]); Figura Soluția problemei cu date inițiale asociată acestui sistem va fi considerată > with(DEtools) : > DEplot(sysl Eql,sysl Eq ,xl(t),x (t),t= , > [[xl( )=l,x ( )=l]],xl= ,x = , ,scene= [xl(t),x (t)]); x (t) &-Г b xl(t) Figura Calculul simbolic al soluțiilor sistemelor de ecuații liniare > with(DEtools): > DEplot d(sysl Eql,sysl Eq ,xl(t),x (t),t= , > [[xl( )=l,x ( )=l]],xl= ,x = , ,scene= [t,xl(t),x (t)]); Figura Se observă că, funcția de plotare plot afișează curbele plane aq = rri(t) și x = în același sistem de coordonate Funcția with(DEtools') : DEplot, odată cu rezolvarea sistemului afișează perechile de puncte (rri(t), rr (t)) care corespund domeniului de variație a variabilei independente t Funcția with(DEtools') : DEplot d permite reprezentarea in trei dimensiuni a curbei spațiale ce reprezintă soluția sistemului considerat Sistemul de trei ecuații diferențiale liniare cu coeficienți constanți omogen: tfi = — aq —izq X' = — X' —Хз Хз = -Хз ( ) Pentru acest sistem considerăm condițiile inițiale aq( ) = , tr (O) = , хз ( ) = , determinăm soluția problemei cu date inițiale și apoi o vom reprezenta grafic: CAPITOLUL > sys Eql:=diff(xl(t), t)=-xl(t)-x (t) ; sys -Eql := ^-fxl (t) = — xl (t) — x (t) > sys Eq :=diff(x (t) ,t)=-x (t)-x (t) ; sys -Eq := ^-tx (t) = — x (t) — x (t) > sys Eq :=diff(x (t) ,t)=-x (t); sys Eq := -^хЗ (t) = —x (t) > dsolve({sys Eql,sys Eq ,sys Eq },{xl(t),x (t),x (t)}); { xl (i) = / (JC — C t + Cl)e~t, x (t) = — ( C — C ) e~l, x (i) = C в-* } > dsolve({sys Eql,sys Eq ,sys Eq ,xl( )= ,x ( )= , x ( )= }); { xl (i) = / ( i + )e-* x (t) = — te t, x (t) = e t} > plot([l/ *( *t~ + )*exp(-t),- *t*exp(-t), *exp(-t)], t=- ,colour=[green,black,blue],thickness=[ , , ], style=[line,point,line]); Figura Sistemul de două ecuații diferențiale liniare cu coeficienți constanți neo-mogen: f = ^i- ^ + t fiisl x- = — rri— x + + t Calculul simbolic al soluțiilor sistemelor de ecuații liniare Pentru acest sistem considerăm condițiile inițiale Xi(l) = , x (l) = , determinăm soluția problemei cu date inițiale reprezentăm grafic acaestă soluție: > sys Eql:=diff(xl(t),t)=xl(t)-x (t)+ *t~ ; sys Eql := -jjxl ( ) = xl ( ) — x ( ) + > sys Eq :=diff(x (t),t)=- *xl(t)- *x (t)+ + *t; sys Eq := ^x ( ) = — x ( ) — x ( ) + + > dsolve({sys Eql,sys Eq }); { xl ( ) = e~ t C + e t Cl - x (t) = ^e~ t-C - e t Cl + + , } > dsolve({sys Eql,sys Eq ,xl( ) = ,x (l)= }); { xl ( ) = у e~ e t — / e e~ t — t , x ( ) = — у e~ e t — / e e~ t + + t } > xl:= / *exp(- )*exp( *t)- / *exp( )*exp(- *t)-t~ : > x :=- / *exp(- )*exp( *t)- / *exp( )*exp(- *t)+ *t+ *t~ : > plot([xl,x ],t=- ,color=[red,green],style= [line,point]); Figura Capitolul Teoreme de existentă și unicitate Metode numerice Proprietăți calitative ale soluțiilor Integrale prime Teoreme de existentă si unicitate pentru ecuatii diferentiale de ordinul întâi neliniare Fie problema cu date inițiale X = f (t,x); x(to) = xo ( ) cu f : Q C IR > IR și (t , x ) G Q Vom enunta și demonstra o teorema referitoare la existența unei soluții (locale) a problemei cu date initiale ( ) Consideram in acest scop doua numere a > si b > , astfel ca dreptunghiul A = { (t, x) | \t — t si t e Ih avem (t,xn(t)) e Q Folosim metoda inductiei matematice, vom arata ca pentru orice t e Ih avem (t,xn(t)) e A Etapa I (a verificării): Pentru n = avem: Xi(t) = Xo + / f(t,Xo(t))dT; to si deci |x (t) — xo| si (V)t G Ih avem (t,xn(t)) G A Trecem acum sa evaluam maximul modulului |x„+ (t) — x„(t)| pe Ih Pentru aceasta, tinem seama de egalitatile: x„+i(t) = xo + / f(t,x„(t))dT to x„(t) = xo + f (t,x„ i(t ))dT to pe care le scadem si obținem: |xra+i(t) — x„(t)| De aici rezulta ca: |x(t) — y(t)| Pentru t G J, t fixat, trecem la limita pentru e si obtinem ca x(t) = y(t), (v)t G J Problema Se stie ca materia radioactivă se dezintegreaza si viteza de dezintegrare este Teoreme de existență si unicitate pentru ecuații diferențiale de ordinul întâi proporțională, la orice moment, cu cantitatea de materie radioactiva rămasă Dacă x(t) reprezinta cantitatea de materie radioactiva ramasa atunci x = —a • x, unde a este o constanta pozitivă în cazul dezintegrarii carbonului radioactiv C , a = si deci în acest ’ caz, ecuatia devine x = -— • x Aratati ca, daca la un moment to se cunoaste cantitatea de carbon radioactiv C dintr-o mostra de animal sau planta gasită într-un strat geologic, atunci se poate reconstitui varsta acelei mostre Rezolvare Fie xo cantitatea de carbon radioactiv C dintr-o mostra la momentul to Problema cu date initiale: x= x(to) = xo •x are solutie unica si aceasta este dată de x(t; t ,x ) = xoe- (t-t ) Daca x este valoarea normala a cantitatii de carbon radioactiv C în starea vie a animalului sau plantei atunci, egaland x = xoe вою(t t ) gasim o ecuatie în t, care ne da timpul t în care animalul sau planta erau vii, iar diferenta t — to arata varsta mostrei Problema Un rezervor cilindric are o gaura circulara la baza prin care lichidul din rezervor se poate scurge O întrebare asemanatoare cu cea din Problema este urmatoarea: daca la un moment dat vedem ca rezervorul este gol putem oare sa știm daca acesta a fost odată plin si cand? Raspunsul este evident nu Cum se explica? CAPITOLUL Rezolvare: Fie x(t) înalțimea lichidului din rezervor la momentul t Notam cu A aria bazei cilindrului si cu a aria gaurii După legea lui Toricelli avem: a ,— r- x = —T\/ q • Cx A Daca I este înalfimea rezervorului, atunci x = I corespunde la situatia cand rezervorul este plin si x = la situatia cand rezervorul este gol Daca la momentul t = rezervorul este plin, atunci x( ) = I si avem: a ,— • t wuu = de unde: I Timpul de golire este t* = — / ^ si dech a \ g x(t) = A \ pentru t* rezervorului daca acesta a fost plin la reprezintăa legea de golire momentul t = Exista o infinitate de solutii x(t) ale ecuatiei a Г x = —тх/ q •Ox A care pentru t = t* sunt egale cu zero Acestea sunt date de formula: ( q • a —-țț- (t* — t — т) pentru A XT (t) = — т t* — т * a Solutia xT (t) reprezinta legea de golire a rezervorului care la momentul — т a fost plin Pe lângă aceste solutii problema Cauchy x = A x x(t*) = Teoreme de existență si unicitate pentru ecuații diferențiale de ordinul întâi mai are ca soluție funcția identic nula x(t) = Aceasta soluție corespunde situației in care rezervorul nu a fost umplut niciodata Rezulta astfel ca solutiile problemei Cauchy descriu toate situatiile posibile și multitudinea acestora se manifesta prin neunicitatea solutiei problemei Cauchy Concluzii Daca variabila de stare x(t) a unui proces fizic sau chimic este solutia unei probleme cu date initiale X = f (t,x), x(t ) = x , atunci aceasta variabila de stare trebuie căutată printre solutiile acestei probleme Cauchy Daca problema cu date initiale in cauza are o singură solutie, atunci gasind-o este clar ca aceasta este cea care descrie evolutia in timp a variabilei de stare Daca problema cu date initiale în cauza are mai multe solutii, atunci gasind una din aceste solutii nu avem nici un drept sa sustinem ca aceasta este aceea care descrie evolutia în timp a variabilei de stare Mai precis, avem nevoie de informatii suplimentare care să permita identificarea acelei solutii care descrie evolutia variabilei de stare în caz de neunicitate gasirea unei solutii a problemei cu date initiale nu înseamna rezolvarea problemei de fizica CAPITOLUL Teoreme de existentă si unicitate pentru sisteme de ecuatii diferentiale de ordinul întâi neliniare Definiția Un sistem de ecuații diferențiale neliniare de ordinul întâi explicit este un sistem de n relații de dependență funcțională de forma: x = f (t,x ,x , ,Xn) x = f (t,Xi,X , ,Xn) ( ) Xn = fn(t, X ,X , ,Xn) dintre un sistem de n funcții necunoscute x , x ,xn și derivatele acestora x ; x ,x,, n în sistemul ( ) funcțiile f , f ,fn sunt funcții reale considerate cunoscute definite pe I x D; I C IR , I interval deschis si D C !Rn, D domeniu Definitia Un sistem ordonat de n funcții reale x , x , , xndefinite pe un interval J C I, de clasa C este o soluție a sistemului ( ) dacă (t, x (t), , xn(t)) E I x Q, (V)t E J si: ' dx dt f (t, x (t), x (t), , xn(t)) f (t, x (t), x (t), , xn(t)) dxn = fn(t, x (t), x (t), , xn(t)) dt pentru orice t E J Definitia Fiind date: t E I si (x , , xn) E D problema determinării soluției x (t), , xn(t) a sistemului ( ) care verifica xi(t ) = x pentru i = ,n, se numeste problemă cu date inițiale sau problemă Cauchy Pentru reprezentarea matriceala a sistemului ( ) notam F : I x D IRn functia matriceala definita prin F(t,x ,x , ,xn) = (f (t,x , ,xn), f (t,x , ,xn), ,fn(t,x , ^n)) Teoreme de existență si unicitate pentru sisteme de ecuații diferențiale de ordinul întâi și cu X matricea coloană (xi,x , ,xn)T Cu aceste notații sistemul ( ) se scrie sub forma matriceala: X = F (t,X) ( ) In aceasta problema derivarea functiei matriceale X(t) inseamna derivarea elementelor matricei Observația Problema cu date initiale (problema Cauchy) se scrie matriceal sub forma: f X = F(t,X) ( ) [ X(to) = X si consta în determinarea functiei matriceale X(t) care verifica ( ) si condifia initialaă X (t ) = X Observația O functie X : J С I R de clasa C este solutie a problemei ( ) daca X = F(t,X(t)), (V)t G J si X(t ) = X (se presupune ca t G J) Consideram a > , b > astfel ca cilindrul A: A = {(t, x) : |t — t ; M = max IIF (t, X )|| si ț M K + J ( ,ж)еД K este constanta lui Lipschitz: ||F(t,X') — F(t,X")|| , de unde se obtine: \\X'(t) - X"(t)|| Trecand la limita pentru se obtine egalitatea X'(t) = X"(t), (V)t G J, t — fixat Proprietăți calitative ale soluțiilor Proprietăți calitative ale solutiilor Consideram I C IR un interval deschis, D C Rn un domeniu, F : I x D ]Rn o functie vectoriala si sistemul de ecuatii diferențiale de ordinul întâi scris sub forma matriceala: X = F (t,X) ( ) Fie X : J C I D si X : J C I D doua solutii locale ale sistemului ( ) Definiția Zicem că soluția locală X este o prelungire a soluției locale X si notăm X facandu-se analog Consideram familia de functii {ха}аеЛ, xa : Ia IR , t G Ia, formata cu toate solutiile locale ale problemei: {X = f (t, x) ( ) x(t ) = X Teorema Cauchy-Lipschitz de existenta si unicitate a solutiei locale pentru problema cu date initiale în cazul unei ecuatii diferentiale de ordinul întai, asigura faptul ca aceasta familie de functii are cel putin un element CAPITOLUL Consideram intervalul deschis I = |J Ia, precum și funcția x = x(t) definita pe I în modul următor: ”pentru t G I consideram a G Л, astfel ca t G Ia si definim x(t) = xa(t)” Aceasta definitie este corecta daca, pentru orice t G Ia П Ia", avem xa' (t) = xa" (t) Analizam aceasta implicare pentru t > to (cazul t t , ti G Ia П Ia", astfel ca xa (t ) = xa" (t ) Consideram in continuare t* = inf{ti : ti > to,ti G Ia' C Ia" ,xa (ti) = xa" (ti)} si remarcam ca xa' (t* ) xa" (t*) = x * Punctul (t*,x*) este în domeniul Q, (t*,x*) G Q, si putem considera constantele pozitive a*, b*, K* astfel ca: - dreptunghiul Д* = {(t,x): |t — t*| astfel ca e Notam cu I = (а ,в ) intervalul de definitie a solutiei saturate a problemei cu date initiale: x = f (t,x), x(t ) = x(t ; t ,x ) = x ( ) Deoarece x(t ; t ,x ) = x , functia x = x(t; t ,x ) este solutie locala a problemei cu date initiale ( ) Rezulta ca I С I si x(t; t , x ) = x(t; t , x ), (V) t G I Pe de alta parte, din faptul ca x(t; t , x ) este solutie saturata, rezultă ca I D I și x(t; t , x ) = x(t; t , x ) Teorema Dacă sunt îndeplinite urmatoarele conditii: (i) e - ); (ii) sirul de numere {tn}n din I converge la e (respectiv a ); Proprietăți calitative ale soluțiilor (iii) șirul de vectori {X(tn; t ,X )}n este convergent la un vector Л, atunci punctul (в , Л) (respectiv (a , Л)) aparține frontierei domeniului Q = I x D; (в , Л) G dQ (respectiv (a , Л) G dQ ) Demonstrație: Facem demonstrația pentru n = urmând ca pentru n > sa se faca în mod analog Pentru orice t G I punctul (t,x(t; t ,x )) aparține domeniului Q si, prin urmare, punctul (в , A) apartine aderentei domeniului Q, (в , A) G Q Putem arata ca (в , A) apartine frontierei dQ, aratand că (в , A) nu apartine domeniului Q Rationam prin reducere la absurd si presupunem ca (в , A) G Q Consideram a > , b > , K > , M > astfel ca dreptunghiul Д = {(t, x) : |t — в | astfel ca e si ti t + e > в , rezulta ca solutia saturata x(t; t , x ) este prelungibilaă, ceea ce este absurd Teorema Dacă I = IR , D = !Rn și в —то), atunci solutia saturata X(t; t ,X ) este nemărginita pe [t ,^ ) (respectiv (a , t ])- CAPITOLUL |x(t; t ,x ) - x , astfel ca |x(t; t , x )| se rationeaza in mod analog Teorema Dacă I = IR , D = IRn și F este lipschitziană pe orice banda de forma Д = J x IRn, unde J C IR este un interval compact oarecare, atunci orice soluție saturată este definită pe IR Demonstrație: Pentru cazul n = , rafionăm prin reducere la absurd si presupunem ca intervalul de definire I = (а ,в ) al solutiei saturate x = x(t; t ,x ) este marginit la dreapta: в Consecința Daca I = (a, b), D = IRn si в a), atunci solutia saturata este nemărginită pe intervalul [t , в) (respectiv (O ,t ]) Consecința Daca I = (a, b), D = IRn F este lipschitziana în raport cu X pe orice bandă de forma J x I, unde J C IR este un interval compact inclus în I, atunci orice solutie saturată este definita pe I Proprietăți calitative ale soluțiilor Teorema Daca funcția de clasa C , F : IR x !î'n IRn nu depinde de t, atunci pentru orice X G R și t , t G IR , avem: X(t + t ; , X ) = X(t ; , X(t ; , X )) = X(t ; , X(t ; , X )) Demonstrație: Demonstram teorema pentru n = Se observa ca pentru t = au loc: x(t + t ; , x ) = x(t ; , x(t ; , x )) = x(t ; , x(t ; , x )) în continuare se remarca faptul ca avem egalitatea: d —x(t + t ; , X ) = f (x(t + t ; , X )) dt si deducem ca x(t +t ; , x ) este solutia saturată a problemei cu date inițiale X = f (x), x( ) = x(t ; ,x ) Rezulta în acest fel egalitatea: x(t + t ; , x ) = x(t; , x(t : , x ), pentru orice t în particular, pentru t = t , se obtine prima egalitate din enunt Pentru cazul n > teorema se demonstreaza analog Fie I un interval real deschis (I G -R ), Q un domeniu deschis în IRn (Q G Rn) si F : I x Q Rin, F = F(t, X) o functie de clasa C Consideram în continuare conditia initiala (t , X ) G I x Q si solutia maximala X = X(t; t , X ) a problemei Cauchy X = F (t,X), X (t )= X Notam cu I intervalul de definire al solutiei maximale X = X(t; t ,X ) Teorema (de dependența continua de ”poziția” inițiala X ) Pentru orice interval compact I* = [T ,T ], inclus în intervalul I (I* C I ), O care contine punctul t în interior (t GI*) și pentru orice e > , exista CAPITOLUL = astfel ca, pentru orice X cu ||X — X || și bt > , astfel ca cilindrul At = {(т, X) : |t — t| , astfel ca d(Y, Z) > r, pentru orice Y G Г si Z G d( j Atj) j= Tubul de securitate A, definit prin: A = {(t,X): t G I* si ||X — X(t; t ,X )|| , astfel ca, pentru orice (t,X ), (t,X ) G A sa avem: ||F(t,X ) — F(t, X )|| De aici rezulta ca, cel putin pentru unul din numerele a , a , definite prin ai =inf {tGI*nIi : ||X(t; t ,X ) —X(t; t , X )|| T Din nou raționăm prin reducere la absurd și admitem de exemplu ca в , există = (e,I*), astfel ca, pentru orice conditie initiala (t ,X ) cu \t - to\ și K > , astfel ca tubul Д, definit prin Д = {(t,X) : t G I*, ||X - X(t; t°,X°)|| | Rezulta de aici ca cel putin pentru unul din numerele a , a definite prin a = inf {t G I* CI : ||X (t ; t°,X ) - X (t ; t°,X°)^ T Rationam prin reducere la absurd si presupunem de exemplu ca в , există S = S(e,I*) > , astfel încat, daca ||p — p || si r > astfel ca multimea Д si S definite prin: Д = {(t,X) : t G I*, ||X — X (t; t ,X ,p )|| astfel încat sa avem: ||F(t,X ,p) — F(t, X , p)|| Rezulta de aici ca, cel putin pentru unul dintre numerele a , a definite prin: Г E ai = inf G I*ПIM : ||X(t; t ,X ,p) — X(t; t ,X ,p )\\ T Raționam prin reducere la absurd presupunand de exemplu в există S = S(e,I*) > astfel încat dacă \\p — p \\ Y + • ek Funcția Hk este continua in raport cu t G Ig, lipschitziana în raport cu Y pe intervale compacte I С I In plus, functia Hk(t,t ,X ,p ,h,Y) tinde la Hk(t, t , X , p , , Y) pentru h uniform pe compacte în raport cu (t, Y) Problemele Cauchy: dYk dt = Hk(t,to,X ,p ,Yk) Yhk (to) = dY k dt = Hk(t,to,X ,p ,Yk) Yk(to) = au solutii definite pe Ig si lim Yfi (t) = Yk(t) uniform în raport cu t pe orice h^ interval J С I Pe de alta parte se verifica usor ca pentru h = avem: ykp\ tvpt vo , ih ) vp + v” ,, ) yh (t) — — |x (t; to, X , p + h • e ) — X (t; to, X , p )] h si deducem ca functia X(t; t ,X , p) are derivate partiale în raport cu pk în (t,to,Xo,p ) si d dt (t,to,X ,p )} \opk J Hk д X t,t ,X ,p , ,— (t,to,X ,p ) dpk (to,to,X ,p ) dpk CAPITOLUL Pe de alta parte: lim Hk(t,t ,X ,p , ,Y) = Hk(t,t ,X°,p°, ,Y) p ^p° uniform în raport cu (t,Y) pe mulțimi compacte dX dX Deci -—(t, t , X , u ) tinde la -— (t, t , X , u ) pentru и u° uniform дрк дик pe orice interval compact I С I Aceasta demonstraza ca derivatele partiale in raport cu pk ale functiei X (t; t , X , u) sunt continue în raport cu p în (t; t , X , p ) în plus, avem: d dt (^X (t,t ,X ,p ^ = \dpk ) д F (t X (t-t V() /, ) ,, \ dX (t-t V() // i dF (t V(t-t V() // ) /, дХF (t, (t; t , , u у г\ (t; t , , u у I r\ (t, (t; t , ,u ),u ) дрк дрк x — (t j t ,X°,u°) = дрк Prin urmare functia X = X (t; t , X , p) este diferentiabila în raport cu p în punctul (t,t ,X ,p ) si pentru orice t G I avem: d dt(дpX(t; t ,X ,p ^ = дХF (t, X (t; t , X , u ), u ) ’ дpX (t; t ,X )+ дpF (t,X (t; t ,X ),p ) дpX(t ; t ,X ,p ) = Metode numerice Metode numerice Metoda liniilor poligonale a lui Euler de determinare numerică locala a unei solutii neprelun-gibile în cazul sistemelor diferentiale de ordinul А A • întai Fie I un interval real deschis și nevid I C IR , D un domeniu nevid în spațiul IRn, D C R si F o functie de clasa C , F : I x D IRn Pentru (t ,X ) G I x D consideram soluția neprelungibila X = X(t; t ,X ) a problemei cu date inițiale X = F (t, X); X (to) = X ( ) si notam cu Io = (ао,в ) C I intervalul de definitie al acestei solutii O prima metoda de determinare numerica locala a solutiei neprelungibile X = X(t; t , X ) este cea a liniilor poligonale a lui Euler In cele ce urmeaza vom prezenta aceasta metoda Consideram doua constante pozitive a si b, a > ,b > astfel ca cilindrul Д = {(t,X)| |t - t | si alegem b = b(e) astfel ca pentru orice (t,X'), (t",X") G Д care verifica inegalitatile: |t' — t"| , alegerea unui J(e) > cu proprietatea menționata este posibila Consideram acum un numar pozitiv h > pe care-l vom numi pas si pe care îl alegem astfel încat sa satisfaca inegalitatea to + a Pentru i = , q consideram numerele ti = t + ih si t-i = t — ih Aceste numere definesc o diviziune a segmentului Ia t — a este egal mărginită pe Ia = [t — a, t + o;] Egalitatea: X£(t) =X° + J F(r, X£(r')')dr + J £(r)dr, (V) t G Ia to to împreună cu inegalitatea: \\Ѳ£(т)\\ h:= : n:= : > f:=(t,x)->-x(t)+ *exp(t): > t:=(n,h)->n*h: > x:=proc(n,h); > if n= then x( ) else x(n- ,h)+h*f(t(n- ,h),x(n- ,h)) end if; > end proc: > x( ):= : > x(t):=[seq(x(i,h),i= n)]; x (t) := [ , , , , , , , , , , ] > t:=[seq(t(i,h),i= n)]; t := [ , , , , , , , , , , ] Metode numerice Pentru a compara valorile numerice ale soluției date de procedura de iteratie Euler cu cele obtinute prin calcul simbolic, vom calcula solutia (obtinuta in Capitolul ) în diferite puncte: > sol x(t):=exp(t)+exp(-t): > eval(sol x(t),t= ); eval(sol x(t),t= ); eval(sol x(t),t= ); eval(sol x(t),t= );eval(sol x(t),t= ); Se observa ca rezultatele obtinute prin calcul numeric cu procedura de iteratie Euler sunt apropiate de cele obtinute prin calcul simbolic doar pentru valori ale lui t apropiate de conditia initiala (zero) aceasta intrucat domeniul de convergenta este mic In concluzie, metoda liniilor poligonale a lui Euler ne da o buna aproximare a solutiei doar pe intervale mici In continuare, vom prezenta un alt exemplu în care vom determina numeric (utilizand procedura de iteratie Euler) solutia sistemului de ecuatii diferentiale: {aii = —xi+ x ( ) X = X + X solutie care a fost deja determinata prin calcul simbolic în Capitolul : xi(t) := ' e t 'e X (t) := - t t ' e t + • e Programand în Maple se obtine: CAPITOLUL > h:= : n:= : > f :=(t,x ,x )->-x (t)+ *x (t): f :=(t,x ,x )->x (t)+x (t): > t:=(n,h)->n*h: > x :=proc(n,h); > if n= then x ( ) else x (n- ,h)+h*f (t(n- ,h),x (n- ,h),x (n- ,h))end if; > end proc: > x :=proc(n,h) > if n= then x ( )else x (n- ,h)+h*f (t(n- ,h),x (n- ,h),x (n- ,h))end if; > end proc: > x ( ):= : x ( ):= : > x (t):=[seq(x (i,h),i= n)]; xl (t) := [ , , , , , , , , , , ] > x (t):=[seq(x (i,h),i= n)]; x (t) := [ , , , , , , , , , , ] > t:=[seq(t(i,h),i= n)]; t := [ , , , , , , , , , , ] Pentru a compara valorile numerice ale soluției cu cele obținute prin calcul simbolic, vom calcula valorile soluției obținute în Capitolul în cateva puncte: > sol x := / *exp( *t)- / *exp(- *t): > sol x := / *exp( *t)+ / *exp(- *t): > eval(sol x (t),t= ); eval(sol x (t),t= );eval(sol x (t),t= ); Metode numerice eval(sol x (t),t= );eval(sol x (t),t= ); > eval(sol x (t),t= ); eval(sol x (t),t= ); eval(sol x (t),t= ); eval(sol x (t),t= ); eval(sol x (t),t= ); Și în acest caz se observa ca, rezultatele obținute numeric cu procedura de iteratie Euler sunt apropiate de cele obtinute prin calcul simbolic doar pe un interval mic Deci, si în cazul sistemelor de ecuatii diferentiale, metoda liniilor poligonale a lui Euler ne da o buna aproximare a solutiei doar pe intervale mici Metoda Runge-Kutta de determinare numerică a unei solutii neprelungibile în cazul sistemelor diferențiale de ordinul întai Cea mai raspandită metoda de determinare numerica a unei solutii neprelungibile este metoda lui Runge-Kutta Ea a fost pusa la punct la sfarsitul sec al XlX-lea de matematicienii germani C Runge si W Kutta In esenta, este tot o aproximare a solutiei neprelungibile cu linii poligonale Convergenta catre solutia saturata este însa mult mai rapida decat în cazul liniilor poligonale a lui Euler Aceasta datorita modului de alegere a ”pantei” In practica, sunt folosite cateva forme particulare ale acestei metode: metoda Runge-Kutta de ordinul al doilea rk , al treilea rk , al patrulea (standard) CAPITOLUL rk și repectiv cea de ordinul al cincilea Fehlberg-Runge-Kutta rkf Fara a intra in detalii privitoare la convergenta, redam aici procedura de iterație a metodei Runge-Kutta standard rk : ti+i = ti + h Xi+ = Xi + h • mR-K unde mR-K = -(mi + m + тз + mA mi = F (ti,Xi) m = F(ti + h/ ,Xi + h/ • m ) m = F(ti + h/ ,Xi + h/ • m ) m = F(ti + h, Xi + h • m ) Aceasta procedura de trecere de la (ti,Xi) la (ti+ ,Xi+ ) este usor de programat Pentru exemplificare vom considera aceleasi exemple ca si in cazul metodei Euler, dupa care vom compara rezultatele numerice obtinute Programand în Maple procedura de iteratie Runge-Kutta standard rk corespunzatoare ecuatiei diferentiale ( ) obtinem: > h:= : n:= : > f:=(t,x)->-x(t)+ *exp(t): > t:=(n,h)->n*h: > x:=proc(n,h) local k ,k ,k ,k ; > if n= then x( ) else k :=f(t(n- ,h),x(n- ,h)); k :=f(t(n- ,h)+h/ ,x(n- ,h)+h*k / ); k :=f(t(n- ,h)+h/ ,x(n- ,h)+h*k / ); k :=f(t(n- ,h)+h/ ,x(n- ,h)+h*k ); x(n- ,h)+h/ *(k + *k + *k +k ) Metode numerice > end if; > end proc: > x( ):= : > x(t):=[seq(x(i,h),i= n)]; x (t) := [ , , , , , , , , , , ] > t:=[seq(t(i,h),i= n)]; t := [ , , , , , , , , , , ] > sol x(t):=exp(t)+exp(-t): > eval(sol x(t),t= ); eval(sol x(t),t= ); eval(sol x(t),t= ); eval(sol x(t),t= ); eval(sol x(t),t= ); Comparând aceste rezultate numerice cu cele obținute cu metoda Euler și apoi cu cele obtinute prin calculul simbolic (din paragraful precedent), observam ca metoda lui Runge-Kutta are domeniul de convergenta întreg intervalul considerat, solutiile obtinute prin rk fiind foarte apropiate de cele obfinute prin calcul simbolic Programand în Maple procedura de iteratie Runge-Kutta standard rk corespunzatoare sistemului de ecuatii diferentiale ( ) obtinem: > h:= : n:= : > f :=(t,x ,x )->-x (t)+ *x (t):f :=(t,x ,x )->x (t)+x (t): > t:=(n,h)->n*h: > x :=proc(n,h) local k ,k ,k ,k ; > if n= then x ( ) else k :=f (t(n- ,h),x (n- ,h),x (n- ,h)); k :=f (t(n- ,h)+h/ ,x (n- ,h)+h*k / ,x (n- ,h)+h*k / ); CAPITOLUL k :=f (t(n- ,h)+h/ +h*k / ,x (n- ,h)+h*k / ); k :=f (t(n- ,h)+h/ ,x (n- ,h)+h*k ,x (n- ,h)+h*k ); x (n- ,h)+h/ *(k + *k + *k +k ) > end if; > end proc: > x :=proc(n,h) local m ,m ,m ,m ; > if n= then x ( ) else m :=f (t(n- ,h),x (n- ,h),x (n- ,h)); m :=f (t(n- ,h)+h/ ,x (n- ,h)+h*m / ,x (n- ,h)+h*m / ); m :=f (t(n- ,h)+h/ ,x (n- ,h)+h*m / ,x (n- ,h)+h*m / ); m :=f (t(n- ,h)+h/ ,x (n- ,h)+h*m ,x (n- ,h)+h*m ); x (n- ,h)+h/ *(m + *m + *m +m ) > end if; > end proc: > x ( ):= : x ( ):= : > x (t):=[seq(x (i,h),i= n)]; x (t) := [ , , , , , , , , , , ] > x (t):=[seq(x (i,h),i= n)]; x (t) := [ , , , , , , , , , , ] > t:=[seq(t(i,h),i= n)]; t := [ , , , , , , , , , , ] > sol x := / *exp( *t)- / *exp(- *t): > sol x := / *exp( *t)+ / *exp(- *t): > eval(sol x (t),t= ); eval(sol x (t),t= );eval(sol x (t),t= ); Metode numerice eval(sol x (t),t= );eval(sol x (t),t= ); > eval(sol x (t),t= ); eval(sol x (t),t= ); eval(sol x (t),t= ); eval(sol x (t),t= ); eval(sol x (t),t= ); Și în cazul sistemului considerat se observa o buna convergența a metodei Runge-Kutta rk Calculul numeric al solutiilor unor ecuatii diferențiale si sisteme de ecuatii diferentiale în aceasta sectiune, utilizand metodele de calcul numeric pe care ni le ofera programul Maple — în care sunt incluse programele procedurilor de iteratie Euler sau Runge-Kutta, vom determina solutiile unor ecuatii diferentiale și respectiv sisteme de ecuatii diferentiale, dupa care, vom reprezenta grafic aceste solutii Primul exemplu se refera la ecuatia diferentiala neliniara de ordinul al treilea cu coeficienti variabili: t 'x + tx + x = lnx, x> ; ( ) > eq :=t~ *diff(x(t),t,t,t)+ *t*diff(x(t),t,t)+ *diff(x(t),t) =ln(x(t)) ; eq := t x (t) + ^ x (t) + dx (t) = ln (x (t)) > dsolve({eq ,x( )= ,D(x)( )= / ,(D@@ )(x)( )= }); Deoarece Maple nu afiseaza nimic înseamna ca este incapabil de a gasi o solutie utilizand calculul simbolic; mai precis, nu poate exprima solutia CAPITOLUL problemei cu date inițiale folosind funcții elementare în acest caz, vom rezolva numeric aceasta ecuatie folosind o sintaxa dsolve care sa permita rezolvarea ecuatiei printr-una din metodele numerice clasice: metoda linilor poligonale a lui Euler, metoda Runge-Kutta de ordin doi, trei sau patru, etc Noua sintaxa dsolve/numeric/classical (numerical solution of ordinary differential equations), specifica calculului numeric, are una din urmatoarele forme: dsolve(odesys, numeric, method=classical); dsolve(odesys, numeric, method=classical[choice], vars, options); in care: odesys numeric method = classical vars options - ecuatia sau lista de ecuatii si conditiile initiale - nume care indica lui dsolve sa rezolve prin metode numerice - optionala, se indica numele metodei numerice: impoly pentru metoda liniilor poligonale a lui Euler; rk , rk , rk pentru metoda lui Runge-Kutta de ordin doi, trei, patru,etc - lista de variabile dependente (optionala) - diferite optiuni: output-ul care dorim s,a se afiseze, numarul de puncte, etc Daca nu folosim optiunea in care sa specifiăm o metoda numerica atunci, calculatorul va alege metoda Fehlberg-Runge-Kutta de ordinul cinci (method=rkf ) — metoda cu cea mai rapidă convergenta Ecuatia ( ) se rezolva numeric cu metoda rk astfel: > a:=dsolve({eq ,x( )= ,D(x)( )= / ,(D@@ )(x)( )= },numeric, method=classical[rk ],output=listprocedure): > sol x := subs(a,x(t)): > sol x( );sol x( );sol x( );sol x( ); sol x( );sol x( );sol x( );sol x( ); Metode numerice Prin instrucțiunea dsolve/numeric/classical, Maple a calculat valorile numerice ale solutiei ecuatiei considerate in punctele domeniului de definiție Pentru afișarea valorilor solutiei x(t) în t = , , , , , , , s-a folosit functia subs Pentru reprezentarea grafică a solutiei ecuatiei diferentiale rezolvata numeric se utilizeaza o functie de plotare specifica metodelor numerice, care are urmatoarea sintaxa: odeplot(dsn, vars, range, options); în care: dsn - numele output-ului ecuatiei rezolvata numeric vars - variabila independenta si functia care se ploteaza (optional) range - optional options - numarul de puncte, diferite modalitati de afisare a solutiei Folosind aceasta instructiune de plotare pentru reprezentarea grafica (Figura ) a solutiei ecuatiei ( ) se obtine: > with(plots):odeplot(sol x,t= ,numpoints= ); CAPITOLUL ao Figura Al doilea exemplu este ecuația diferențiala liniara de ordinul al doilea: r d i „ di L • dt + R ' dt + Ci = -E» • ' ( ) a carei soluție i(t) exprima intensitatea curentului într-un circuit R-L-C (vezi Capitolul ) Considerând valori numerice pentru R, L, C, E si reducand ecuatia la un sistem de doua ecuatii diferentiale se obtine: > R:= : C:= : L:= : E= : v:=Pi/ : > Eq:=L*diff(i(t),t,t)+R*diff(i(t),t)+( /C)*i(t)=-E*w *sin(t*w) : > sys Eq :=diff(x (t),t)=x (t): > sys Eq :=diff(x (t),t)=- *x (t)- *x (t)- *Pi/ *sin(t*Pi/ ): Pentru rezolvarea numerica a acestui sistem si pentru vizualizarea solutiilor corespunzatoare la diferite conditii initiale (Figura si Figura ) se folosesc functiile: with(DEtools) : DEplot with(DEtools) : phaseportrait with(DEtools) : DEplot d > with(DEtools):DEplot({sys Eq ,sys Eq },{x (t),x (t)}, t= ,[[x ( )= ,x ( )= ],[x ( )= ,x ( )= ], [x ( )= ,x ( )= ]],scene=[t,x (t)],method=classical[rk ]); Metode numerice xl(t) : - - : Figura > with(DEtools):DEplot({sys Eql,sys Eq },{xl(t),x (t)}, t= ,[[xl( )= ,x ( )= ],[xl( )=l,x ( )=l], [xl( )= ,x ( )= ]],scene=[t,x (t)],method=classical[rk ]); ; Ѳ:- x (t) ~ - ■ : - Figura In aceste figuri sunt reprezentate soluțiile (rri(t), ^г(і)) pentru trei condiții inițiale Se observă că, indiferent de condițiile inițiale, după un anumit timp soluția se stabilizează in jurul unei soluții periodice Acest fapt, reiese și din portretele de fază care se obțin cu: > with(DEtools):phaseportrait([sys Eql, sys Eq ], [xl(t),x (t)],t= ,[[xl( )= ,x ( )= ]], scene=[xl(t),x (t)],method=classical[rk ]); CAPITOLUL Figura > with(DEtools):phaseportrait([sys Eql,sys Eq ], [xl(t),x (t)],t= ,[[xl( )=l,x ( )=l]], scene=[xl(t),x (t)],method=classical[rk ]); Figura Portretele de fază (Figura și Figura ) arată faptul că, soluțiile sistemului (intesitatea curentului și variația acesteia) se stabilizează după un anumit timp, tinzănd către un ciclu limită Mai precis, folosind condiții inițiale din interiorul sau exteriorul ciclului limită soluțiile se stabilizează in jurul unei soluții periodice Același fenomen de stabilizare se observă și din figura tri-dimensională (Figura ): > with(DEtools):DEplot d({sys Eql,sys Eq }, {xl(t),x (t)},t= ,[[xl( )= ,x ( )= ]], scene=[t,xl(t),x (t)],method=classical[rk ]); Metode numerice Figura Un alt exemplu este sistemul lui Lotka-Volterra de două ecuații diferențiale neliniare care constituie un model matematic utilizat in biologie care descrie evoluția m timp a două specii pradă-prădător (de exemplu sardine-rechini): ( ) unde x(t) reprezintă numărul sardinelor, iar z/(i) numărul de rechini Cu instrucțiunile with(DEtools') : DEplot (Figura și Figura ) și respectiv with(DEtools') : DEplot d (Figura ) se obține evoluția m timp a celor două specii: > with(DEtools):DEplot({diff(x(t),t)=x(t)*(l-y(t)), diff(y (t),t) = *y(t)*(x(t)-l)},{x(t),y(t)}, t= ,[[x( )= ,у( )= ]],scene=[t,x(t)], linecolor=t/ ,method=rkf ); CAPITOLUL : : : к© ; : : Figura > with(DEtools):DEplot({diff(x(t),t)=x(t)*( -y(t)), diff(y(t),t)= *y(t)*(x(t)-l)},{x(t),y(t)}, t= ,[[x( )= ,y( )= ]],scene=[t,y(t)], linecolor=t/ ,method=rkf ); ; : ; У© Ș : о в; ; Figura > with(DEtools):DEplot d({diff(x(t),t)=x(t)*(l-y(t)), diff(y(t),t)= *y(t)*(x(t)-l)},{x(t),y(t)},t= , [[x( )= ,y( )= ]],scene=[t,x(t),y(t)] , stepsize= ,linecolor=t/ ,method=rkf ); Metode numerice Figura Portretul de fază a evoluției celor două specii este prezentată in Figurile : > with(DEtools):DEplot([diff(x(t),t)=x(t)*( -y(t)), diff(y(t),t) = *y(t)*(x(t)-l)],[x(t),y(t)] ,t= , [[x( )=l ,y( )=l ],[x( )=l,y( )= ],[x( )= ,y( )= ]], stepsize= ,title='Lotka-Volterra model', color=[ *y(t)*(x(t)-l) ,x(t)*(l-y(t)) , ] , linecolor=t/ ,arrows=MEDIUM,method=rkf ); CAPITOLUL Lotka-Volterra model : : : : b : : Figura Acest sistem are soluțiile staționare ( , ), ( , ) și soluții periodice (Figura ) Interpretarea acestora este următoarea: i) soluția staționară ( , ) reprezintă dispariția ambelor specii; ii) soluția staționară ( , ) reprezintă situația de echilibru (numărul de sardine este egal cu numărul de rechini); iii) soluțiile periodice care înconjoară soluția staționară ( , ) reprezintă variații ale numărului de sardine și respectiv rechini intre două limite Aceste variații (creșteri sau descreșteri) descriu lipsa sau abundența de hrană (sardine) care duce la micșorarea sau creșterea numărului de prădători (rechini) Un alt portret de fază interesant este al sistemului: {i = у — z у = z-x ( - ) c = x — y, Metode numerice care arată că, din orice punct ar pleca soluțiile, toate vor evolua către zero (Figura ): > with(DEtools):phaseportrait([D(x)(t)=y(t)-z(t), D(y)(t)=z(t)-x(t),D(z)(t)=x(t)-y(t)* ],[x(t),y(t),z(t)], t=- ,[[x( )= ,y( )= ,z( )= ]],stepsize= O , scene=[z(t),y(t)],linecolour=sin(t*Pi/ ), method=classical[rk ]); y(t) r-i—i—|— —i ’i' i' i—i—rr - - - - - Figura CAPITOLUL Integrale prime Consideram sistemul de ecuații diferențiale neliniare de ordinul întâi explicit: X = f (t,Xi,X , ( ) Xn fn (t, X , x , •••, xn) în care funcțiile f , f ,fn sunt funcții reale de clasa C definite pe I x D; I C IR , I - interval deschis si D C IRn - domeniu Definiția Se numește integrală primă a sistemului ( ) o funcție U : I x D IR de clasa C care nu este constantă, dar este constantă pe soluțiile sistemului ( ) Teorema Condiția necesară si suficienta pentru ca o functie U : I x D IR de clasă C neconstanta să fie integrală primă este ca U să verifice: n + dT ’ fi = (V) (t,x , ,xn) G I x D Demonstrație: Fie (to,x , ,хП) G IxD si x (t),xn(t) solutia sistemului ( ) care verifica x^t ) = x , i = ,n si U(t, x (t),xn(t)) o integrala prima Deoarece U(t, x (t),xn(t)) este constanta, rezulta ca dU, , — (t,x (t), ,xn(t)) = dt De aici avem ca dU n dU (t,x (t), ,x„(t)) +J ^ (t,x (t), ,xn(t))• fi(t,x (t), ,xn(t)) = , (V) t G I Pentru t = , rezulta de aici egalitatea: dU (t xo xo), dU (t xo xo) f (t xo xo) dt (to, x , •••, xn) + / y dx (to, x , •••, xn) fî(to, x , •••, xn) = / г Integrale prime întrucât (t , x , ,хП) este un punct oarecare din mulțimea I x D rezulta dU n dU I V f = (V) (t x x ) G I x D dt f г , (V)(t,x ,•••, xn) G x D• i= Sa aratăm acum implicatia reciproca, adica: dacă U: IxD IR este o functie de clasa C , care nu este constanta si verifica dU dU ~dt + ~dt • fi = atunci U este constanta pe solutiile sistemului ( ) In acest scop consideram o solutie oarecare x (t), ,xn(t) a sistemului ( ) si functia ^(t) = U(t,x (t), ,xn(t)) Pentru a arata ca functia ,xn, c , ,cm)j ,рт(( xm+ , c , ,cm) , xm+ , ,xn) dt dxn , problema naturala care se pune este: care este numărul maxim de integrale prime independente pentru sistemul ( )? Raspunsul la aceasta întrebare este data de urmatoarea teorema Teorema Pentru orice punct (to,x , x^,) G IxD exista o vecinătate deschisă V С I x D astfel ca pe V sistemul ( ) are n integrale prime independente și orice integrală prima definită pe V se exprimă ca funcție de cele n integrale prime independente Demonstrație: Fie X(t; t ,X ) solutia saturata a sistemului ( ) care verifica X(t ) = X si Io intervalul ei de definitie Consideram un interval compact [ , T ] inclus în Io care contine punctul t în interior Conform teoremei de continuitate în raport cu conditiile initiale exista o vecinatate Uo a lui x = (x , , x^) astfel încat pentru orice X' = (x/, , xn') G Uo solutia sistemului ( ) care coincide cu X 'în t = t este definita pe [ , ]; notam cu X(t; t ,X') aceasta solutie saturata Functia (t, X') Л ■ X(t; t ,X') este de clasa C pe baza teoremei de diferentiabilitate în raport cu conditiile initiale Din aceeasi teorema rezulta d ca matricea ——[X(t; t , X')] este nesingulara si prin urmare aplicatia X' ■ dX' X(t; t , X') este local inversabila Exista deci doua vecinatati deschise W , W ale punctului (t ,X ) si o functie fi de clasa C , fi : W prietatile: fi(t,X(t; t ,X')) = (t, X') si X(t; t ,fi(t,X")) = X" W si (t,X") G W Componentele scalare fik ale functiei fi raman constante x , ,xn se înlocuieste cu o solutie a sistemului definita în vecinatatea considerata, deci fik sunt integrale prime In plus fi(t ,X') = (t ,X') de unde rezulta ca rang () = n Adica integrabile prime fi , ,fin sunt inde- \dx; / pendente W cu pro- (V) (t,X') G atunci când CAPITOLUL Fie acum U o integrala prima oarecare a sistemului ( ) Conform definiției U(t,X(t; t ,X')) nu depinde de t si putem scrie U(t,X(t; to,X')) = h(X') înlocuind X’ cu ^(t,X") si ținand seama de egalitatea X(t; to,^(t,X")) = X" obtinem: U (t,X") = h(^(t,X")) Aceasta din urma egalitate arata ca integrala prima U este o functie h de cele n integrale prime independente ^ , ,^n- O metoda de determinare a integralelor prime este data de ecuatia: f + t • fi = în aceasta ecuatie U este functie necunoscuta si intervine in ecuatie prin intermediul derivatelor partale de ordinul întâi De aceea ecuatia aceasta se numeste ecuatie cu derivate partale de ordinul intai Orice solutie a acestei ecuatii este o integrala prima Observam ca daca functiile fi ,fi , ,fin sunt astfel ca n fio + fii • fi = i dU dU si exista o functie U cu proprietatea — = fi și —— = fii, atunci functia U dt dx* este integrala prima pentru sistemul ( ) Integrale prime Exerciții: Fie sistemul de ecuații: X = x I X = —x Sa se determine o integrala prima [ Цо = , Ц = X , //- = X R: U (X ,X ) = (x + x ) în descrierea miscarii solidului rigid intervine sistemul: A • p = (B — C)g • r B • q = (C — A)r • p C • Г = (A — B)p • q Sa se determine doua integrale prime pentru sistemul considerat R: U (p,q,r) = Ap + Bq + cr si U (p,q,r) = A p + B q + C r Sa se determine doua integrale prime pentru sistemul: dX dy X y dz —z R:U (x, y, z) = X^/y si U (x, y, z) = xz b) dX dy dz z y X z y X R: U (x,y, z) = x + y + z si U (x, y, z) = x + y + z c) dx dy dz x (y + z) —y (z | x) z (y — x) R: U (x, y, z) = xyz si U (x, y, z) = —|- — xyz CAPITOLUL d) dx dy dz xy x y z(x + y ) R: Ui(x,y,z) = x - y și U (x,y,z) = X■ Capitolul Ecuații cu derivate parțiale de ordinul întai Ecuații cu derivate parțiale de ordinul întai liniare Definitia O ecuație cu derivate parțiale de ordinul întâi liniara este o relație de dependența funcțională de forma: n ~ EUU , , — • fi(X ,X , ,Xn) = • n dXi i= ( ) înțre derivațele parțiale de ordinul înțâi ale funcției necunoscuțe u = u(x , x- , , xn) Funcțiile f , f , , fn în ecuația ( ) fi : Q C Rin+ IR se consideră cunoscuțe si sunț presupuse funcții de clasa C pe Q Definitia O soluție a ecuației ( ) esțe o funcție u : Q C Q IR de clasa C asțfel încaț în țoațe puncțele (x , x , , xn) G Q să avem verifi-cațăa idențițațea: n du ^ (xo,x , ,xn) • fi(x ,x , ,xn) = i= i CAPITOLUL Fie C = (C , £i, , Cn) € Q astfel încât n ^f (Co,Ci, -,Cn) = o i= Existenta unui punct C € Q cu aceasta proprietate poate fi admisa pentru câ, in caz contrar, toate funcțiile fi ar fi identic nule pe Q de unde ar rezulta ca ecuatia ( ) este verificata, oricare ar fi functia u de clasa C pe Q Putem admite de asemenea ca f (C ,C , ■ ■ ■ ,Cn) = Aceasta intrucat din n ipoteza fi(C , C , ■ ■ ■, Cn) = rezulta ca exista i € { , , , n} astfel ca i= fi (C) = Ceea ce presupunem noi în plus este faptul ca i = , adica f (C) = - Daca f (C) = , atunci se face rafionamentul pentru fi Revenim deci la f (C) = si deoarece f este continua, exista o vecinatate V a punctului C astfel ca f (x) = pentru orice x € V• Consideram acum sistemul de ecuatii diferentiale: xi = gi(t,x , ■ ■ ■ ,xn), i = ,n ( ) unde /, fi(t,x , ■ ■ ■ ,xn) gi(t, x , ■ ■ ■ , xn) ț ,, x f (t,x , ■ ■ ■ ,xn) t fiind componenta x a vectorului x = (x , x , ■ ■ ■, xn) € V, t = x Sistemul ( ) este definit pentru (t,x , ■ ■ ■ ,xn) € V,iar functiile gi sunt de clasa C pe Vț Teorema Soluțiile ecuației ( ) definite pe V C Vț sunt integrale prime ale sistemului ( ) și reciproc Demonstrație: Fie u : V R o solutie a ecuatiei ( ) Cu notatiile introduse avem: dU XC/ X dt(t,x , ■ ■ ■ ,xn) • f (t,x , ■ ■ ■ ,xn) + n n EdU , \ , X —(t,x , ■ ■ ■ ,xn) • fk (t,x , ■ ■ ■ ,xn) = dxk k= k Ecuatii cu derivate parțiale de ordinul întâi liniare deci (t x ) + V — (t x, x ) • h= dt (t,x„ ,x„)+ dxt (t,x!, -,xn) /o(t,x , ,x„) - și prin urmare: du , \—л du , , \ d- (t,xi, ,Xn) + (t, xi, ,Xn) • gfc(t,xi, ,Xn) = k=i xk ceea ce arata ca u este integrala prima a sistemului ( ) Cu un rationament asemanator se arata ca daca u este integrala prima pentru sistemul ( ) atunci este soluție pentru ecuația ( ) Prin urmare, problema determinarii solutiilor ecuatiei ( ) revine la problema determinaarii integralelor prime pentru sistemul ( ) Ținand seama de cele demonstrate pentru integralele prime ale unui sistem, rezulta ca solutiile ecuatiei cu derivate partiale ( ) au urmatoarele pro-prietaati: oo o oo o i) oricare ar rr (xo, xi, , xn) g □ □, daca jo(xo, xi, , x n) , atunci exista o vecinatate deschisa V C Q a punctului (# , ,#£ ) si n functii ui, ,un : V R de clasa C astfel încât: a) ui, u , , un sunt solutii ale ecuatiei ( ); b) uk(x ,xi, ,xn) = xk; c) det f = , k,l = , , ,n ii) daca u este o solutie oarecare a ecuatiei ( ) definita pe V C V, atunci exista o vecinatate deschisa D a lui (я , ,яП) si o functie Y : D C R Ri astfel ca u(xo, xi, , xn) = y(ui (x , xi, , xn) , u (x , xi, , xn) , , un (x , xi, , xn) ) Definitia (Problema Cauchy pentru ecuația ( )) Se numește problema Cauchy pentru ecuația ( ) problema determinării unei soluții u a ecuatiei ( ) astfel ca u(x , xi, , xn) = h(xi, , xn), (V)(xi, , xn) G D' C D, unde x° si functia h : D C IRn Ri de clasă Ci sunt date; D fiind o vecinătate deschisă a lui (x , ,x)() CAPITOLUL Teorema Daca în (x , x°, , x^) G Q avem /o(x ,x , ,x°) = = atunci pentru h de clasa C definită într-o veci problema Cauchy are soluție unică nătate deschisă D a lui (xi, , x ) Demonstrație: Consideram solutiile u , u , în vecinatatea deschisa V a lui (x , xl, , x^ ,un ale ecuației ( ) definite ) cu proprietatea ufc(x ,x , ,xra) = xk, k = , ,n Exista o vecinătate V С V a lui (x , x , , astfel ca pentru (x , x , , xn) G V sa avem (u (x , x , , xn), , un(x , x ■ , , xn)) G D- Fie u(x , x , , xn) = h(u (x , x , , xn) , , un(x , x , , xn)) Functia u(x , x , , xn) definita astfel este ( ): de clasa C si verifica ecuatia r\ n П r\ r\ du r du dh oui x + dxk k dyi x dh f dul dyi у x n n л r\ dh dul f ăyi'd^k' k'- n \ f + E g-A = k= k / în plus, u(x , x , , xn) = h(u (x , x , , xn) = h(x , ,xn) , ,un(x ,x , ,xn)) = deci u este solutia problemei Cauchy Pentru unicitate sa presupunem ca u este o a lta solutie a problemei Cauchy definita pe V Rezulta ca exista astfel ca u(x , x , , xn) (ui (x , xi , • • • , xn) , • • • , un (x , xi , • • • , xn) )• Ecuatii cu derivate parțiale de ordinul întâi liniare De aici rezulta ca h(xi, ,Xn) = u(x , X , ,Xn) = y(xi, ,Xn); și deci u coincide cu u Daca în punctul (x , x , , x^) avem fk (х ,х , ,хП) = , atunci se poate formula un rezultat analog pentru o functie h definita pe o T ПО fof Ci ЛІ CXQ /~»Ті ICO Q ТЛ ITIZ^fillllll ( dXn Pentru A = rezulta: X -(X ,X , ,Xn) + + Xn^-(X ,X , , Xn) = , (V)(X ,X , , Xn), dX dXn de unde avem ca: / х х u(X ,X , ,Xn) = Y I , \ xn xn xn— Xn CAPITOLUL Exerciții: Rezolvați următoarele ecuații cu derivate parțiale de ordinul întâi liniare: a) b) du du y • — x • p- = dx dy R: u(x,y) = y (x + y ) du du x • + у • л- = dx dy R: u(x,y) = Y c) X • -z— y ' дУ - z ' d~z = R: u(x,y,z) = Y (xy/y,xz) d) du n -г ( du du\ du xy-t — v - У\ У • Г~ z • ^ = xy • dx \ dy dz) dz R: u(x, y, z) = Y \ yz + xu y + ^/ — y \ , x • earcsiny Rezolvați urmatoarele probleme Cauchy: du du a) x • ——+ y • — = ; u(x, ) = x dx dy x R: u(x, y ) = Ecuatii cu derivate parțiale de ordinul întâi liniare du du du b) Vх • + Ѵу + ^z ■ = ; u(l,y,z) = y - z dx dy dz R: u(x,y,z)= y + z — y/yz c) ( + x ) ’ дс + ХУ ’ dj = ; u( , У') = У R: ( ) y u(x'y) = T+X- CAPITOLUL Ecuatii cu derivate parțiale de ordinul întâi cvasiliniare Definiția O ecuație cu derivate parțiale de ordinul întai cvasiliniara este o relație de dependență funcționala de forma: n - • fi(x , ,xn,u) - g(x , ,Xn,u) = i ( - ) dințre funcția necunoscuța u = u(x , ,xn) și derivațele parțiale ?—,i = , n ale acesțeia Funcțiile f , , fn și g în ecuația ( ) se considera date: fi,g :Q C R ' R si sunt presupuse de clasa C Ecuația ( ) se numeste cvasiliniara pentru ca este liniara in derivatele parțiale ——, dar nu este liniara in general in u Definiția O soluție a ecuației ( ) este o funcție: u : D C Rn —> R de clasa C care are propriețațea ca pențru orice (x , ,xn) G D puncțul (x , , xn, u(x , , xn)) G Q si n du (x , ,xn) • fi(x , ,xn,u(x , ,xn))- g(x , ,xn,u(x , ,xn)') = în D Pentru determinarea solutiilor ecuatiei ( ) consideram ecuatia cu derivate partiale de ordinul întai liniara: n dv dv — • fi(x , ,xn,u) + — • g(x , ,xn,u) = ( - ) cu (x , , xn, u) G Q Ecuatii cu derivate parțiale de ordinul întâi cvasiliniare Teorema Dacă v este soluție a ecuației ( ) defintă pe Q' C Q, iar (X , , хП, u ) G Q' și dV(x°, ,Xn,u ') = atunci funcția u = u(xi, , xn) definita implicit prin ecuația v(x , , xn, u) — v(x°, , хП, u ) = este o solutie a ecuatiei ( ) Demonstrație: Fie u = u(x , ,xn) funcția definita implicit prin ecuația v(x , , xn, u) — v(x , , хП, u ) = Avem: v(x , , Xn, u(x , , Xn)) — v(x , , xn, u ) = și prin derivare în raport cu xk gasim: dv -—(xi, ,Xn,u(xi, ,Xn)) + dv du — (xi, ,Xn,u(Xi, ,Xn)) • (Xi, ,Xn) = du dxk pentru k = , , , n înmulțind pe rând aceste egalitati cu fk(xi, , xn, u(x , , xn)) si adunandu- le gasim: n Я EăX"(X’“(X)) • fk(X,u(X))+ k=l k A(X,u(X)) (X,u(X)) • fk(X,u(X)) = u k=i xk unde X = (x , , xn) Dar v fiind solutie a ecuatiei ( ) are loc n V Țț—(xi, ,Xn,u(xi, ,Xn)) • fk(xi, ,Xn,u(xi, ,Xn)) + y dXk k=i k CAPITOLUL dv —(xi, ,xn,u(xi, ,xn)) • g(xi, ,xn,u(xi, ,xn)) = și prin urmare: dv — (xi, ,xn,u(xi, ,xn)) • g(xi, ,xn,u(xi, ,xn)) = xi, •••, xn) ' fk (xi, •••, xn, u(xi, xn) ) • n —(xi, ,xn,u(xi, ,xn)) k=i dv Tinem seama de faptul ca — (xi, , xn, u(xi, , xn)) = deducem egali-du tatea: n du У'' -—(xi, ,xn) • fi(xi, ,xn,u(xi, ,xn)) = g(xi, ,xn,u(xi, ,xn)) y dxk k=i k Rezulta in acest fel ca functia u = u(xi, ,xn) este solutie a ecuatiei ( ) Teorema arata ca rezolvarea ecuatiilor cvasiliniare cu derivate parțiale de ordinul întâi se reduce la rezolvarea ecuatiilor cu derivate partiale de ordinul întâi liniare Definiția Pentru (x , , xn) G Q problema găsirii unei soluții u(xi, , xn) a ecuației ( ) astfel încat u(x ,x , ,xn) = ((x- , ,xn), ț fiind o funcție data de clasa Ci se numeste problemă Cauchy pentru ecuația ( ) Teorema Daca fi(x , , xn, ((x , , x^) = , atunci există o vecinătate deschisa V a punctului (x , ,x°fi) și o solutie u = u(xi, ,xn) a ecuatiei ( ) definită pe V astfel încat u(x ,x , ,xn)= //'? / R: y (x + y , y (u + Vu + )) = Rezolvați următoarele probleme Cauchy: du du y - x a dx dy u u( ,y) = y b) R: u (x, y, z) = (x + y - ) + xy - (x + y - ) du du xy -y • — = x, u( ,y) = — dx dy y R: u(x, y) = x + xy Ecuatii cu derivate parțiale de ordinul întâi cvasiliniare c) u • - — + (u — x ) • ^ + x = , u(x, x ) = x dx dy R: y + u (x, y) = (x • u(x, y) — y) CAPITOLUL Calculul simbolic al soluțiilor ecuațiilor cu derivate parțiale de ordinul întai Pentru calculul simbolic al soluțiilor ecuațiilor cu derivate parțiale de ordinul întai Maple folosește funcția pdsolve (find Solutions for parțial differential equations (PDEs) and systems of PDEs) cu una din urmatoarele sintaxe : pdsolve(PDE, f, INTEGRATE, build); pdsolve(PDE system, funcs, other options); în care: PDE - ecuația cu derivate parțiale pe care dorim să o rezolvam f - numele funcției necunoscute (implica existența derivatelor parțiale ale mai multor funcții) INTEGRATE -(opțional) indica integrarea automata a ecuațiilor diferențiale ordinare care intervin atunci cand PDE este rezolvatăa cu metoda separăarii variabilelor build - (opțional) indica afisarea unei forme explicite (daca este posibil) a soluției Pentru exemplificare, consideram ecuația cu derivate parțiale de ordinul întai liniaraă: du du У x = dx dy ( ) Cu instrucțiunea pdsolve se obține soluția generala (toate soluțiile) a ecuației: > PDE := y*diff(u(x,y),x)-x*diff(u(x,y),y) = ; PDE := ydu (x, y) — xdu (x, y) = > pdsolve(PDEl); u (x, y) = -F (x + y ) Calculul simbolic al soluțiilor ecuațiilor cu derivate parțiale de ordinul întâi Se observa ca, soluția generala este afișată cu ajutorul unei funcții F care poate fi orice funcție de clasa C Pentru determinarea unei anumite soluții avem nevoie de o conditie initiala, dar in sintaxa functiei pdsolve prezentata anterior nu exista nici un parametru care sa specifice utilizarea acesteia în cele ce urmeaza, vom mai determina solutia generala pentru o ecuatie cu derivate partiale de ordinul întai în care functia necunoscuta este de trei variabile x, y, z: du у 'ir + дУ ( ) > PDE :=sqrt(x)*diff(u(x,y,z),x)+sqrt(y)*diff(u(x,y,z),y)+ sqrt(z)*diff(u(x,y,z),z) = ; PDE := y/xdu (x, y, z) + ^/yЦu (x, y, z) + sqrtzu (x, y, z) = > pdsolve(PDE ); u (x, y,z) = x^-x- - /у} Capitolul Ecuații cu derivate parțiale de ordinul al doilea liniare Clasificarea ecuațiilor cu derivate parțiale de ordinul al doilea liniare (clasificarea problemelor de fizică - matematica) Definiția O ecuație cu derivate parțiale de ordinul al doilea liniară este o relație de dependența funcțională de forma: ( ) dințre funcția necunoscuță u = u(xi, ,xn), derivațele parțiale de ordinul ânțâi și de ordinul al doilea ale acesțeia Funcțiile a-ij, bi,c, f : Q C lRn IR din ecuația ( ) se considera cunoscuțe și sunț presupuse cel puțin conținue pe Q Definiția O soluție clasică a ecuației ( ) esțe o funcție u : D C Q IR de clasa C care are propriețațea ca pențru orice (x , , xn) G D avem: n n Clasificarea ecuațiilor cu derivate parțiale de ordinul al doilea liniare +c(xi, ,Xn) • u(xi, ,xra) + f (xi, ,xra) = Fie T : Q Qz un difeomorfism de clasa C de componente țk = țk(x , , xn), k = , n si u = u(x , , xn) o soluție a ecuației ( ) Consideram functia v = u o T- si din derivarea functiilor compuse obtinem: n d u n — dv înlocuind derivatele dțk dți dv dțk dți dxi dxj dțk în ecuatia ( ) obtinem: d țk nn d v n Et dv bk • dT + C V + f = ( ) unde: i= j= nn S & b bi dxi + ai dxfixj i= i= j= Ecuatia ( ) este echivalenta cu ecuatia ( ), solutiile u ale ecuatiei ( ) obtinandu-se din solutiile v ale ecuatiei ( ) cu formula u = v o T Pe de alta parte intr-un punct oarecare, fixat (x) , ,xn) G Q putem considera forma păatraticaă n n • Уі • yj i= j= unde aj = aij(x) , , хП) Prin schimbarea de variabile forma patratica considerata se transforma in forma pătratică: nn nn I dțk dțl > > I > > aij • • nn k= = \i= j = \ A d v n n c ț av • v CAPITOLUL Prin urmare, coeficienții derivatelor parțiale de ordinul al doilea se schimba la fel cu coeficientii formei patratice sub actiunea transformarii liniare: k p • Пк i Se stie ca, prin alegerea unei transformări liniare adecvate, forma patratica se aduce la forma canonică (adica matricea poate fi adusa la forma diagonala: |а^| = sau si, = daca i = j) Aceasta inseamna că alegand in mod adecvat transformarea fk = fk(x , ,xn), k = ,n, coeficientii derivatelor d v partiale — în punctul fk = fk(x , ,x”), k = , n vor fi egali cu + , — sau d v , iar coeficientii derivatelor partiale o vor fi nuli d ) i=i / ( ) Tipul ecuatiei ( ) în (x , , x)) se determina prin aducerea formei patratice: n n yy i= j= la forma canonica Ecuatia poate fi adusa la una din formele standard prezentate alegand transformarea ёк = ёк(x , ,ёп), k = , n astfel încat transformarea liniara: Уі = Ц(x ’ ’xn) • Пк к= i sa aduca forma pătratica la forma canonica Sa examinam in continuare problema posibilitatii aducerii ecuatiei ( ) la ”forma canonică” (una din formele prezentate) într-o vecinatate a unui punct (x , ,xn), in conditiile in care în toate punctele acestei vecinatati ecuatia apartine aceluiasi tip Pentru a aduce ecuatia ( ) la forma canonica într-un anumit domeniu ar trebui, în primul rând, sa impunem functiilor ёк(x , xn), k = , n conditiile diferentiale akl = pentru k = l Numarul acestor conditii este n^n — si este mai mic decat n (n reprezinta numarul functiilor ёк, k = ,n) daca n ecuatia ( ) nu poate fi adusa la forma CAPITOLUL canonica in vecinătatea unui punct (x- ,, xn) Pentru n = elementele nediagonale ar putea fi anulate in general, insa elementele diagonale ar putea fi diferite de ± , Deci, nici in acest caz, ecuația ( ) nu poate fi adusa la forma canonica în vecinatatea unui punct Doar pentru n = , putem anula unicul coeficient nediagonal si satisface conditia de egalitate a celor doi coeficienti diagonali Aceasta înseamna ca doar in cazul n = putem aduce ecuatia cu derivate partiale la forma canonica pe o vecinatate Exerciții Aduceti la forma canonica ecuatia: â u d u d u du du r a • ~~ + a ■ — + a ' + b • ; + b • ;—+ c • u + J (x,y) = dx dxdy dy dx dy unde aij, bi, c sunt constante R: Scriind ecuatia caracteristica asociata: a a + a = se obtine discriminantul: Д = (a ) — a • a - daca Д > (cazul ecuatiei hiperbolice), atunci ecuatia caracteristica are doua radacini reale: {dy a + V (a ) — a • a dx a dy a — (a ) — a • a dx a Se rezolva aceste ecuatii diferentiale ordinare si se obtine f (x,y) = c , respectiv g(x,y) = c Schimbarea de variabile care ne va duce la forma canonica va fi: f ) in cadranele II si IV; - ecuatia este de tip eliptic (Д R: Ecuatia este de tip hiperbolic (Д > ); forma canonica este: d v dv dv dțdy dț dy ’ c) d u d u ■ + •У 'Sy = ’x’y> R: Ecuatia este de tip eliptic (Д = ni, i = , ; VP si VQ sunt funcțiile vectoriale (gradientii funcțiilor P si Q) definite prin: dP dP VP — —— • ei + —— • VQ = dQ , dQ + - Pentru demonstratie se calculeaza membrul stang al egalitatii ( ) tinand seama de formula de integrare prin parti ( ) si se obtine: [[ P • AQdxidx = !! + dr(^\\dxidx JJq JJq Ldxi \dxij dx \ x J [ P dQd ff dP dQd d + / P • cos ai •—— ds — —— dxidx + ddQ dxi JJn dxi dxi + f P dQd [[ dP dQ d JdQ dx JJa dx dx \dQ dQ —— • cos ai + —— • cos a dxi x ds — dP dP dQ\ , — • — dxidx + —— • —— dxidx = Q RQ , dQ —— • cos ai + —— • cos a dxi dx ds— VP • VQdxidx = Formulele lui Green si formule de reprezentare în două dimensiuni [ P •^Qds — [[ VP •VQdx dx dQ dn JJQ Teorema (cea de-a doua formulă a lui Green) Daca Q este un domeniu mărginit în IR cu frontiera dQ netedă (parțial neteda), iar funcțiile P,Q : Q IR sunt de clasă C pe Q si de clasa C în Q atunci are loc următoarea egalitate: ([ (PXQ — QAP)dxidx =[ (p • — Q • ^ds ( ) Jq J dQ \ dn dn / Demonstrație: Egalitatea ( ) se obține scriind egalitatea ( ) pentru funcțiile P • VQ si Q • VP dupa care se face diferenta acestora Teorema (de reprezentare a unei functii de două variabile) Daca Q este un domeniu marginit în R cu frontiera dQ netedă (parțial netedă) §i functia u : Q IR este de clasa C pe Q §i de clasa C în Q atunci pentru orice X = (x ,x ) G Q are loc următoarea egalitate: u(X) X и • Xu(Y )dy dy + I' du + П Lln X Я • -(Y ' — г d f \ — — / u(Y) • -— ( In ————r I dsY L dnY\ ||X — Y\\J Y ( ) Demonstrație: Functia E(X ) = — ln definita pentru orice X G R, X = este de clasa C pe IR \ { } si are proprietatea XE = Pentru X G IR , X-fixat consideram functia: E(X — Y) ln n \\X — Y || definita pentru orice Y G IR \ {X} Aceasta functie de Y este de clasa C si are proprietatea XYE(X — Y) = De asemenea, pentru orice Y G IR , CAPITOLUL Y-fixat functia E(X — Y) = — ln ————— definita pentru orice X G n ||X — Y || IR , X = Y este de clasa C si verifica AXE(X — Y) = Consideram X G Q, X-fixat si e > astfel ca, pentru orice Y cu ||X—Y || rezulta: lim -[ ln^——— •-^-^uțY )dsY = e-o nJdn ||X - YII BUY Y = -[ ln^— • ă—u(Y}dsY n Jdn ||X - Yу } Y CAPITOLUL lim d Л Ъ дЦln dSY u(Y) • n JdQ d f \ ln |jx-Yj[) d + u(X>■ Rezulta de aici ca avem: • (Au)(Y )dyidy = f , du = - X ln X - • -(y + г d f \ — / u(Y) • -— ( ln^———- | dsY + u(X), ^/dQ } dnY\ ||X - Y\\J Y V >’ sau u(X) ln Q цх Yц • (Au)(Y + • ' JKr |X -Й (Y)dsY- -[ u(Y) • —— ( ln - ) dsY nJdQ^ dnY\ ||X - Y\\J Y Comentariu: Aceasta formula de reprezentare permite determinarea funcției u în toate punctele lui Q cunoscand urmatoarele elemente: Laplacianul du functiei u pe Q, valorile derivatei normale — pe frontiera dQ si valorile dn functiei u pe frontiera dQ Vom arata in continuare importanta acestei formule de reprezentare pentru cazul funcfiilor armonice Formulele lui Green si formule de reprezentare în două dimensiuni Definitia Zicem că o funcție u : Q R este armonică în Q dacă funcția u este de clasă C si Au = , (V)(x , x ) G Q unde Ли reprezintă Laplacianul funcției u: d u d u u x + dx ' Observația în condițiile din teorema de reprezentare ( ) dacă funcția u este armonica in Q și de clasa C pe Q atunci avem: u(X) = ! І^л, (Y)dsy - n IIх — Y|| dny г d f i \ — / u(Y) • -— ( In ————r | dsy - ny\ ||X — Y\\J Y Teorema (de reprezentare a funcțiilor armonice) Daca u : Q C Ri Ri este o funcție armonică pe domeniul Q atunci (V)X G Q si (V)r > pentru care discul: B(X,r) = {Y : ||Y — X|| astfel încat B(Xo,ro) = {Y : ||Y - X || astfel încat pentru orice X G B(X ,r ) avem: u(X) Aplicam formula de reprezentare a lui u(X ) pe dB(X ,p) si gasim: u(X ) = -— [ u(Y )dsY np J dB(X ,p) Frontiera dB(X ,p) se descompune astfel: dB(X , p) = dB(X , p) П B(X , n) U B(X ,p) П (Q \ B(X ,n)) = U si cu aceasta descompunere formula de reprezentare devine: u(X ) = -ds /* u(Y)dsY + -[ u(Y)dsY astfel încât pentru orice X G P discul închis: B(X, r) = {Y : ||Y — X|| Formulele lui Green si formule de reprezentare în dimensiunea n > în acest paragraf prezentăm câteva rezultate de calcul diferențial și de calcul integral pentru functii de n variabile (n > ), care intervin in rezolvarea ecuatiilor cu derivate parfiale in dimensiunea n Teorema (de legătură între integrala pe un domeniu mărginit și integrala pe frontiera acestuia) Dacă Q C R este un domeniu marginit cu frontiera dQ neteda (parțial netedă) și f , f ,fn sunt n functii fi, f ,fn : Q IR continue pe Q și de clasa C pe Q, atunci are loc egalitatea: cos(n, ei)dS ( ) Demonstrație: Semnificata simbolurilor din formula ( ) este aceeași ca si în cazul n = , iar demonstrata teoremei se face in mod analog Teorema (de integrare prin pă/rți) Daca Q C IRn este un domeniu mărginit cu frontiera dQ netedă (partial neteda) si f,g sunt două funcții f,g : Q IR , continue pe Q si de clasa C pe Q, atunci are loc egalitatea: [ gdx ••• dxn = [ f • g • cos(n,ei)dS — i f •-^ dx ••• dxn ( ) JQ dxi JdQ JQ dXi Demonstrație: Analoaga cu cea din cazul n = Teorema (prima formulă a lui Green) Daca Q C IRn este un domeniu mărginit cu frontiera dQ netedă (partial neteda) si f,g sunt doua functii f,g : Q IR , de clasă C pe Q si de clasă C în Q, atunci are loc egalitatea: Г f dg f / f • Ag dx ••• dxn = f dS — Vf •Vgdx ••• dxn ( ) Jq Joq dn Jq CAPITOLUL Demonstrație: Simbolurile din formula ( ) au următoarele semnificații: n — og д A fg — - dg —cos(n e-) W = V fe- dn Qxi ' г ' Qxi Teorema se demonstreaza ca si în cazul n = Teorema (a doua formulă a lui Green) Dacă Q C IRn este un domeniu mărginit cu frontiera dQ netedă (parțial netedă) și f,g sunt două functii f,g : Q IR , de clasă C pe Q si de clasă C în Q, atunci are loc egalitatea: [(f Лл — ri Af) dr, dr — f ( — dS ( ) j (I g дІ) dxi axn J g j ao ( ) Demonstrație: Teorema se demonstreaza ca si în cazul n = Teorema (de reprezentare a unei functii de n variabile) Daca Q C IR este un domeniu marginii cu frontiera dQ netedă (parțial neteda) si u : Q IR este o functie de clasă C pe Q si de clasa C în Q, atunci are loc următoarea formulă de reprezentare: u(X) (n Q \\X - Y|| n— Xu(Y) dy • • • dyn + (n \\X - Y||n— dny (Y) dSy ( ) (n - * )^n dQ dny IIX - Y|| n— dSy + d unde an reprezinta aria suprafetei bilei B( , ) = {Y =(y , ,yn) G IRn\\\Y|| Comentariu Formula ( ) permite construcția funcției u din valorile du Laplacianului Au al funcției in Q, valorile derivatei normale — a lui u pe dn dQ si valorile lui u pe dQ Vom arata ce devine aceasta formula de reprezentare în cazul funcțiilor armonice Definiția Zicem că funcția u : Q C lRn IR este armonică în Q dacă funcția este de clasă C în Q și dacă AL d u Au = dff = ’ (V)(xi, ,xra) G Q i— i Observația In condicile din Teorema (de reprezentare), daca functia u este armonica in Q, atunci avem: u(X) = (n - )u„ ||X - Y||n- dnY (Y) S ЭО d ( \ dS dny \J|X - Y\\n- J y Teorema (de reprezentare a funcțiilor armonice) Daca Q C lRn este un domeniu și u : Q R este o funcție armonică în Q (Au = ), atunci pentru orice x G Q si orice r > astfel încât bila închisă B( , r) = {Y G RD\ ||Y|| ( ) iar funcția v(X,Y) are următoarele proprietăți: a) Y v(X, Y) este funcție armonică în Q și continua pe Q pentru orice X = (x , , xn) G Q, X fixat b) pentru X G Q, X fixat, G(X,Y) = , (V)Y = (y , ,yn) G dQ Propoziția Dacă există funcția Green G(X,Y) pentru Problema Dirichlet, atunci ea este unicăa Demonstrație: Daca G , G sunt doua functii Green pentru Problema Dirichlet, atunci pentru orice X G Q fixat, v(X, Y) = V (X, Y) - v (X, Y) este functie armonica in Q si identic nula pe dQ Rezulta v(X, Y) = , (V)Y G Q, si (V)X G Q, fixat Aceasta arata ca v(X, Y) = , adica V (X,Y ) = V (X,Y) Functia Green pentru Problema Dirichlet Propoziția Dacă pentru orice X G Q funcția Y v(X,Y) este de clasă C pe Q, atunci funcția Green este simetrică, adică: G(X ,X ) = G(X ,X ), (V)X ,X G Q Și X = X ( ) Demonstrație: Se considera funcțiile P(Y)=G(X^Y) și Q(Y)=G(X Y) și se aplica cea de-a doua formula a lui Green pentru aceste funcții pe domeniul Q£ = Q \ (B(X ,e) U B(X ,e)), unde B(Xi,e) = {X : ||X-Xi|| Propoziția Daca pentru orice X G Q functia Y v(X,Y) este de clasa C pe Q, atunci functia Green pentru Problema Neumann este simetrică: G(X ,X ) = G(X ,X ), (V)X ,X G Q si X = X ( ) CAPITOLUL Demonstrație: Analoaga cu cea din cazul funcției Green pentru Problema Dirichlet Teorema Daca G este o funcție Green pentru Problema Neumann ( ) și u este o soluție a acestei probleme, atunci există o constantă C astfel ca să aibe loc egalitatea: u(X) = -[ G(X,Y) • f (Y) dY E(X,Y) • g(Y) dSy + C ( ) JQ JdQ Demostrație: Analoaga cu cea de la cazul soluției Problemei Dirichlet Observația Determinarea funcției Green pentru Problema Neumann revine la determinarea functiei v(X,Y) ceea ce înseamna, conform definitiei functiei Green pentru Problema Neumann, rezolvarea Problemei Neumann -Ду v(X,Y) = f, (V)(xi, ,xn) G Q dv dny dE dQ dnY dQ ( ) Aceasta problema este aparent mai simpla decat Problema Neumann ( ), dar în realitate este complexa De aceea demonstrata existentei solutiei Problemei Neumann ( ) prin folosirea functiei Green se poate face doar pentru cazuri pentru care se știe ca există functia Green Probleme pentru ecuatia lui Laplace pe disc Problemele Dirichlet si Neumann pentru ecuația lui Laplace pe disc Separarea variabilelor Consideram mulțimea Q = {(xi, x ) G R x + x G Q u |dQ= h ( ) Solutia acestei Probleme Dirichlet se poate determina folosind functia Green pentru Problema Dirichlet Problema Neumann pentru ecuatia lui Laplace pe discul de raza R înseamna determinarea acelor solutii ale ecuatiei ( ) ale caror derivata dupa directia versorului normalei exterioare la dQ coincide cu o functie continua g dinainte data pe dQ, adica: d u d u d u ,, dxf + ' = , V(xi,x ) G Q dx ( ) =g dQ CAPITOLUL Pentru ca problema ( ) să aibe soluție este necesar ca / g(Y)dsY = ( ) JdQ si daca aceasta condiție este îndeplinită, atunci, folosind functia Green pentru Problema Neumann, se pot determina solutiile prolemei ( ), (care difera printr-o constanta aditiva) Scopul nostru in acest paragraf este sa prezentam o altă metoda, numita ”separarea variabilelor”, cu care se pot determina solutiile problemelor ( ) si ( ) Deoarece s-a considerat Q ca fiind un domeniu circular vom face mai întai o schimbare de coordonate, mai exact vom scrie problemele ( ) si ( ) în coordonate polare: X = r cos p „ Г- n x z ^ , r > , p G [ , n) ( ) X = r sin p Notam cu T transformarea (r,p) (xi,x ) definita de ( ) Daca u este o functie care satisface ecuatia ( ) atunci notam cu u functia u = u o T Folosind regulile de derivare ale functiilor compuse deducem ca u verifica ecuatia: d u £ d u du dr r dp + r dr care se numeste ecuatia lui Laplace în coordonate polare Metoda separarii variabilelor constă în cautarea unor solutii u(r, p) de forma: u(r,p) = P(r) • Q(p) ( ) adica solutii care sunt produse de functii ce depind fiecare de cate o variabila r, respectiv p Impunand la ( ) sa verifice ( ) rezulta: p" P' Q Г — + rp = ( ) Q Membrul stang în aceasta egalitate depinde doar de r iar membrul drept de p Cum r si p sunt variabile independente rezulta ca fiecare membru este constant Daca notăm cu A aceasta constanta atunci deducem din ( ) egalitatile: r P" + rP' - AP = ( ) Probleme pentru ecuatia lui Laplace pe disc Q" + XQ — ( ) Deoarece funcția u este de clasa C soluția ecuației ( ) trebuie să verifice Q( ) — Q( n) De aici rezulta ca Xn — n , n G IV Pentru Xn — n avem: Qn(^) — An cos n^ + Bn sin nv ( ) în care An și Bn sunt constante arbitrare Ecuatia ( ) este liniara de tip Euler si se rezolva făcandu-se schimbarea de variabila r — e* Pentru X — n solutia generala este: Pn(r) — Cnrn + Dnr n, daca n — , , ( ) Po(r) — A + B ln r, daca n — ( ) Rezulta in acest fel ca ( ) admite urmatoarea familie de solutii: Un(r, ^) — Ao + Bo ln r, rn(An cos n^ + Bn sin n^), r n(A n cos n^ + B n sin n^), n— n — , , ( ) n — , , în care Ao, Bo, An, Bn, A n, B n sunt constante oarecare Observația Orice suma finita de solutii de forma ( ) este solutie pentru ecuatia ( ) Definitia O soluție formală a ecuației lui Laplace în coordonate polare este o ”funcție” w(r, ^) de forma: w(r, ^)—Ao + Bo ln r + ^[(Anrn + A nr n) cos v: + (Bnrn + B ,nr n) sin n^] n= ( ) Denumirea de solutie formala provine de la faptul că nu avem informatii relative la convergenta seriei ( ) Ceea ce este clar este ca, termenii seriei ( ) sunt solutii si ca, daca seria are doar un numar finit de termeni diferiti de zero, atunci suma este o functie care este solutie a ecuatiei lui Laplace CAPITOLUL Pentru determinarea soluției Problemei Dirichlet ( ), mai întâi sriem problema în coordonate polare: d u £ d u du dr r дф + r dr ’ u(r, ф) = u(r, ф + n) astfel încat sa avem | an" | restart; > DirichletInt:=proc(f,R) local a ,a,b; > a :=( /( *Pi))*Int(f,phi=-Pi Pi); > a:=n-> /Pi*Int(f*cos(n*phi),phi=-Pi Pi); > b:=n-> /Pi*Int(f*sin(n*phi),phi=-Pi Pi); > a +add(r~n/R~n*(a(n)*cos(n*phi)+b(n)*sin(n*phi)),n= Order); > RETURN(map(simplify,value(%))); > end: Apelarea acestei proceduri se realizeaza cu instrucțiunea Dirichletlnt(f, R) în care f este functia h din problema ( ), dupa cum se poate observa din urmatoarele exemple: Exemplul : d u d u o x ^ dx ’ u(x ,x ~) = (xi + X ) , x + x f:=(x +x )~ : R:=R: > f:=subs(x =r*cos(phi),x =r*sin(phi),r=R,f); f := (Rcos (ф) + R sin (ф)) > sol :=DirichletInt(f,R); soli := R + r sin ( ф) Exemplul : d u d u x + dx ’ u(x , x ) = x • x , x + x f:=x *x : R:= : > f:=subs(x =r*cos(phi),x =r*sin(phi),r=R,f); f := cos (ф) sin (ф) > sol :=DirichletInt(f,R); sol := / r sin ( ф) Exemplul : d u d u du o dr r дф + r dr ’ r f:=sin(phi)~ : R:= : > sol :=DirichletInt(f,R); sol := / r sin (ф) — / r sin( ф) Exemplul : d u £ d u du dr r дф + r dr , r f:=sin(phi)~ +cos(phi)~ : R:= : > sol :=DirichletInt(f,R); sol := / + / r cos ( ф) Capitolul Soluții generalizate Metode variationale în capitolul anterior ne-am ocupat de rezolvarea Problemei Dirichlet și a Problemei Neumann in cazul ecuatiei lui Poisson Particularitatea acestei ecuatii consta în aceea ca ecuatia eliptica este definita de Laplacianul A în cele ce urmeaza vom considera ecuatii eliptice mai generale, numite ecuații eliptice de tip divergența pentru care vom formula conceptul de Problema Dirichlet Pe langa conceptul de solutie clasica introducem si conceptul de solutie generalizata si prezentam conditii în care solutia generalizata exista si este unica Deasemenea , în acest capitol introducem Problema Cauchy-Dirichlet pentru ecuatii parabolice si hiperbolice, si prezentam conditii în care solutia generalizata a acestor probleme exista Ecuatia eliptică de tip divergență si Problema Dirichlet Ecuatia eliptică de tip divergentă si Problema Dirichlet pentru această ecuatie Consideram un domeniu mărginit Q C Rn cu frontiera dQ neteda (parțial neteda) si functiile reale aij, c si F definite pe Q (i,j = , n) avand urmatoarele proprietati: ) aij sunt de clasa C pe Q si exista p > astfel încat pentru orice X = (x , xn) G Q si orice (î , î , , în) G Rn avem n n aij (X) ■ îi ■ îj — îi si aij (X) = aji(X) i,j= i= ) functia c este continua pe Q si c(X) — , (V)X G Q ) functia F este continua pe Q Definiția Se numește ecuație cu derivate parțiale eliptică de tip divergență o relatie de dependenta functională de forma n я / n я \ (X) • âXH + c(X)• u = F(X) (V) X G Q ( ) i= i \j= j / dintre functia necunoscută u și derivatele sale parpiale de ordinul întai și de ordinul al doilea în ecuatia ( ) functiile aij,c si F se considera cunoscute Definitia O functie reala u de clasa C pe Q care verifica n я f n o \ ~ E (X) • fi(X) + c(X) • u(X) = F(X) (V) X G Q i= i j= j se numește solutie clasică a ecuatiei ( ) Definitia Problema determinării acelor solutii u ale ecuatiei ( ) care sunt continue pe Q si verifică conditia la frontieră u|do = h ( ) se numeste Problema Dirichlet pentru ecuatia ( ) CAPITOLUL In egalitatea ( ), h este o funcție reala continua pe , considerata cunoscuta Problema Dirichlet pentru ecuatia cu derivate partiale eliptica de tip divergenta se noteaza cu n o / n Я \ ăV E “v (X) • gx + c(X) •u — F (X )■ (V) X e i= \j= j/ ( ) , u|dQ — h- Observația Daca exista o functie H : ! D IR de clasa C astfel încât H|dQ — h atunci u este solutie a Problemei Dirichlet ( ) daca si numai daca functia v — u — H este solutia Problemei Dirichlet: n n / n Я \ — E (X) • Я + ') • v — G(X), (V) X € i= \j= j / ( ) v|sn — E “ V (X ) • )+ J/ i= \j= J / si n n d(a • u)\ — + c(X) • (a • u) = a • Au, c) I A(a• u)=~Ysdx (aij(X) • i= г \j= ceea ce demonstreaza liniaritatea lui A iii) Pentru u,v G D calculam / Au • vdX si gasim : Jq ( ( \П' д du\ t / Au • v dX = — / —— I У a,j (X) • —— I • v dX + / c(X) • u • v dX Jq Jq ,= dxi \j= dxj j Jq nn dQ i= j= —— cos(u, e,) • v dS+ Ecuatia eliptică de tip divergență si Problema Dirichlet n n Q i= j= du dv f я— я— dX + / OXj OXi Jq • u • v dX n n Q i= j= du dv [ я—я— dX + / dxj dxi • u • v dX Calculând in continuare / Av • u dX gasim: л П П r\ z» ) • fr'TTdX + / dXj dXi J Q = aji(X) rezulta / Au • vdX = Av • udX Q Q i= j = • u • v dX Deoarece aij (X) Teorema Exista y > astfel încât pentru orice u G D sa avem Au•udX >y u dX ( ) Demonstrație: Calculam Jq Au • udX si obtinem : nn du f Y~dX + c(X)u (X) dX > dxj J Q Q i= j = > du dxi dX + c(X)u (X) dX p n Ramane sa evaluam / \ do ~~ dX pentru u G D Aceasta evaluare con- stitue continutul unei teoreme care poarta numele lui Friedrichs Conform acestei teoreme exista o constanta k > , astfel încât pentru orice u G C (Q) cu u|dQ = avem: p n |u(X)| dX ^ • [ |u(X)| dX + [ c(X)u (X)dX > — [ |u(X)| dX Jq k Jq Jq k Jq Astfel inegalitatea ( ) a fost demonstrata Ramane sa demonstram acum teorema folosita în demonstrata teoremei ( ) denumita inegalitatea lui Friedrichs Teorema (inegalitatea lui Friedrichs) Exista o constanta k > , astfel încât pentru orice u G C (Q) cu u|^Q = să avem: |u(X)| dX Deoarece [ |u'(X')| dX' = [ |u(X)| dX si Jq' Jq q dX q i unde u'(X') = u(TX) si X' = TX, rezulta ca este suficient sa se demonstreze inegalitatea ( ) pe Qz Mai mult, Qz fiind domeniu marginit exista a > astfel ca Qz С Га = {X' | CAPITOLUL Rezulta astfel ca u este minim absolut pentru funcționala Фр Reciproc, daca presupunem ca u G D este punct de minim absolut pentru functionala Фр, atunci pentru orice u G D avem Фр (u) > Фр (uo) Consideram u G D, u-fixat, t G IR si remarcam ca: Фр(uo + tu) > Фр(uo) pentru orice t G IR Prin urmare functia ^(t) = Фр (uo + tu) are un minim in t = Din conditia ^'( ) = rezulta: / (Au — F)udX = (V)u G D Jq Deoarece spatiul vectorial D este dens în spafiu Hilbert L (Q) rezulta ca Au = F Observația Aceasta teorema nu este o teorema de existenta pentru ca nu stabileste faptul ca functionala Фр are un punct de minim absolut Teorema reduce insa problema existentei si unicitatii solutiei clasice a ecuatiei Au = F la problema existentiei si unicitatii punctului de minim al functionalei Фр Teorema (teorema de unicitate) Pentru F G C(Q) există cel mult un u G D astfel încat Au = F Demonstrație: Presupunem ca pentru F G C (Q) exista u , u G D astfel încât Au = F si Au = F De aici rezulta Фр(u ) > Фр(u ) si Фр(u ) > Фр(u ) Prin urmare Фр(u ) = Фр(u ) Tinaml seama de inegalitatea: Фр(u ) — Фр(u ) > A |u — u | dX Jq rezulta in continuare ca / |u — u | dX = de unde rezulta egalitatea Jq u = u Ecuatia eliptică de tip divergență si Problema Dirichlet Observația i) Teorema de unicitate demonstrează ca Problema Dirichlet are cel mult o solutie clasica ii) Existența unui minim absolut în D a funcționalei (u) nu este adevarata Exista exemple care demonstreaza ca în general Problema Dirichlet nu are solutie clasica Pentru a asigura minimul absolut al functionalei Ф/ largim spatiul vectorial D astfel ca sa devina spatiu Hilbert In acest scop, pe D x D definim urmatoarea corespondenta: D x D Э (u, v) ^ A= Au • vdX ( ) Jq Lema Corespondența definită cu ( ) este un produs scalar pe spațiul vectorial D Demonstrație: prin verificare Observația Pentru (u,v) G D are loc egalitatea: A= [ VV агз(X)^u^rdX + [ c(X)u(X)v(X)dX Jq i=i =і dxi Jq Definiția Completatul spațiului vectorial D în norma || • ||a generată de produsul scalar A= Au • vdx Jq se numește spatiul energetic al operatorului A și se notează cu XA Observația Elementele spatului energetic XA sunt elementele spatiului vectorial D la care se mai adauga limite în norma || • ||a de siruri fundamentale un din D Altfel spus, daca u G XA atunci u G D sau exista (un)n cu un G D astfel încat ||un — u||A Un element u G XA are proprietatea: u G L (Q), — G L (Q) si u|dQ = CAPITOLUL Observația Spațiul energetic XA împreuna cu produsul scalar A este un spatiu Hilbert Lema Funcționala Фр : D IR' definită de formula (u) = I Au • udX — I F • udX ( ) JQ JQ în care F G L (Q), se prelungește prin continuitate la o funcțională continuă ФP definita pe spafiul energetic XA Demonstrație: Aratam la început ca funcționala Фр (u) definita cu ( ) este continua în topologia spatiului energetic In acest scop consideram u G D, un sir (un)n, un G D cu proprietatea ca ||un — u||A ^- si apoi sirul Фр (un)= Aun • undX — F • undX QQ Aratam ca acest sir numeric este convergent si limita lui este Фр (u) = Au • udX — F • udX QQ Intr-adevar, avem: Фр (un) = ІЫІа — F • undX Q Фр(u) = \\u\\ A — [ F • udX Q de unde rezulta: |фр(un)— фр(u)l а | llunHA — IHA I + /* |FI • lun— u\dX a Q a llunH A — INIA | + \\F IIl (q) • llun— uIIl (q) Deoarece convergenta ||un — u||A implica convergentele I IHIIA — l|uIIAI si \\un — u^L (q) , rezulta convergenta |ФР(un) — ФР(u)| - Ecuatia eliptică de tip divergență si Problema Dirichlet Daca u GD și u G XA atunci se definește ФР(u) cu Фр (u) n n ^ ^ aij(X) i= j= dX F • udX Jq iar pentru un G D, ||un — u||A se reface același raționament Teorema (de existentă și unicitate a punctului de minim absolut) Pentru orice F G L (Q) prelungita ФF prin continuitate a funcționalei ФР la spatiul energetic XA are un singur punct de minim Demonstrație: Pentru a demonstra existenta punctului de minim consideram functionala liniara si continua u i—> F • udX Jq definita pe spatiul energetic XA Conform cu teoremei lui F Riesz exista o functie uF G XA astfel încat sa avem: A= / F • udX Jq pentru orice u G XA Ramane de aratat ca uF este punctul de minim absolut al functionalei ФF(u) în acest scop calculam ФF(uF) si ФF(u) pentru u G XA Avem: Фf(uf) = A — l (q)— ^upHA — ||up||a = — \\uF||A ФF(u) = A — L (Q) = A + A— A— l (q) = \\u — uF ||A — \\uF ||A = \\u — uF ||A + Фp(uF) Aceasta din urma egalitate Фр (u) = ||u — up HA + Ф p (up) este adevarata pentru orice u G XA si arata ca ФF(u) > Фp(up) Egalitatea arata si faptul ca Фp (u) > Фp (up) pentru orice u = uF ceea ce demonstreaza unicitatea CAPITOLUL Observația Daca punctul de minim absolut uF al funcționalei prelungite apartine lui D C XA atunci AuF = F, si prin urmare uF, este solutie clasica a Problemei Dirichlet Observația Daca punctul de minim absolut uF al functionalei prelungite Ф/ - nu apartine spatiului vectorial D C XA atunci nu-i putem aplica operatorul A ca sa vedem daca este solutie a ecuatiei Au = F In acest caz este naturala întrebarea: Ce reprezintă uF în contextul rezolvării ecuației Au = F ? Vom da un raspuns la aceasta întrebare aratand ca operatorul A are o prelungire A : D(A) L (Q), numita prelungirea Friedrichs, unde D C D(A) C XA si ca uF verifica AuF = F Aceasta inseamna ca functia uF nu mai verifica conditia de solutie clasica : u G D si + c(X) • u = F ci verifica doar conditia : uF G XA si f nn n o o f f / aij (X) •^-dx + / c(X) • uf • vdX = / q fi дх dxi F • vdX, (V) v G XA Teorema (de prelungire a lui Friedrichs) Operatorul G : L (Q) XA definit prin G(F)= uF G XA unde uF verifică A= F (Q), (V) uGXA, are urmatoarele proprietăfi: i) G este liniar și marginit ; ii) G este injectiv ; iii) G este autoadjuct ; iv) G este compact ; Ecuatia eliptică de tip divergență si Problema Dirichlet v) G : ImG C XA L (Q) este o prelungire a lui A Demonstrație: Remarcam la început ca operatorul G este corect definit pentru ca existenta si unicitatea functiei uF = G(F) a fost demonstrată i) Pentru а, в G IR , F , F G L (Q) si u G XA avem: A = L (Q) = L (Q) + L (Q) = а L (Q) ' в L (Q) = L (Q) de unde rezulta egalitatea : G(aF + eF ) — aG(F ) + eG(F ) care aratăa căa operatorul G este liniar Pentru a demonstra maărginirea F G L (Q) si evaluam ||G(F)||A Avem: J -||G(F)||b(fi) L (Q) este constanta din inegalitatea lui Friedrichs si jo > este constanta din conditia de elipticitate Simplificand cu ||G(F)||l (q) în inegalitatea obținuta, rezulta: k \\G(F)||L (Q) A= , (V)u GXA L (q)= , (V)u GXa F = iii) Operatorul G : L (Q) ImG C XA fiind injectiv putem considera operatorul G- : ImG L (Q) si aratam ca este autoadjunct, adica : L (Q)= l (q), (V)u,v G ImG Astfel, l (Q) = A= A= l (Q) = l (Q) Folosim acum acest rezultat pentru a demonstra ca operatorul G este au-toadjunct: L W= L (q= = l (Q)= l (Q) iv) Prin faptul ca operatorul G este compact intelegem ca transforma sfera închisa de raza unu din L (Q) intr-o multime compactă din L (Q) Consideram F G L (Q) cu \\F||l (q) l (Q) A= L (q)= A (V)v G XA Rezultă de aici egalitatea : GAu = u, (V) u G D, care demonstrează apartenenta u G ImG si egalitatea G- u = Au Am arătat în acest fel că operatorul G- : ImG C XA L (Q) este o prelungire a operatorului A Definitia Operatorul G se numește prelungirea Friedrichs a operatorului A și se noteaza cu A Teorema (caracterizarea minimului absolut al funcționalei Ф) Funcția uF G XAeste minimul absolut al functionalei ФF : XA IR , ФF(u) = A — l (q) dacă si numai daca uF este solutia ecuatiei Ău = F, unde A este prelungirea Friedrichs a operatorului A Demonstrație: Rezultatul se obtine imediat din construcfia prelungirii Friedrichs a operatorului A Observația Din cele prezentate rezultă că pentru orice F G L (Q) ecuatia Au = F are o singură solutie în spatiul Hilbert XA Definiția Solutia uF a ecuatiei Ău = F se numeste solutia generali-zatăa a ecuatiei Au = F CAPITOLUL Observația Daca uF E D C XA atunci uF este soluție clasica a Problemei Dirichlet Daca uF E XA nu aparține la D (uF țUD) atunci ea verifica doar: n n ti d d f* f* / (x) •-A-F —— dx + c(x) • uF(x) • v(x)dx = F(x) • v(x)dx Jq i=i dxj dxi Jq Jq pentru orice v E XA în continuare vom descrie o metoda de determinare a solutiei generalizate uF a ecuatiei Au = F (despre care stim ca exista si este unica) Metoda se bazeaza pe determinarea valorilor proprii si vectorilor proprii ai prelungirii Friedrichs A Definiția Un număr A este valoare proprie pentru operatorul A dacă există o funcție u în D(A) (domeniul de definiție al operatorului A), u = , astfel încat să avem: Au = A • u Teorema Pentru operatorul A există un sir infinit de valori proprii L (Q) = fij Demonstrație: Operatorul G : L (Q) L (Q) este liniar autoadjunct si complet continuu Pe baza unei teoreme relative la aceasta clasa de operatori liniari rezulta ca, G admite un sir infinit de valori proprii si un sir infinit de functii proprii Valorile proprii ale lui G le notam cu A A Am Ecuatia eliptică de tip divergență si Problema Dirichlet iar funcțiile proprii cu u , u , •••, um, • Presupunem ca aceste valori proprii sunt aranjate in șir astfel ca sa avem: Ai A > ••• > > > iar functiile proprii sunt alese astfel ca sa avem ij• Ținand seama de egalitatea A = G- rezulta ca A admite sirul de valori proprii |A | pentru orice m Intr-adevar din G • um = • um rezulta ca Am m /vm Pe de alta parte, L (Q) = A> llUm^L (Q) si astfel A > llUm|lL Q > Am Aratam acum ca, ăl admite chiar un sir infinit de valori proprii La început aratam ca multimea de definitie ImG a operatorului A ( formata din elementele G(F) cu F G L (Q) ) este un spatiu vectorial infinit dimensional Pentru aceasta, fie u o functie de clasa Ccu suport compact in Q Consideram functia F = Au si observam ca u este solutia clasica a Problemei Dirichlet Au = F Rezulta ca, u este si solutie generalizata a acestei probleme, ceea ce înseamna ca G(F) = u Am aratat în acest fel ca orice functie de clasă C'x cu suport compact inclus în Q apartine multimii Im(G) si ca urmare spatiul vectorial Im(G) este infinit dimensional CAPITOLUL Sa presupunem acum prin absurd ca, operatorul G are doar un număr finit de valori proprii diferite de zero: A , A , , Xm Ținand seama de faptul ca G este autoadjunct si compact rezultă de aici: m G(F) = uk, (V)F g L (Q) k= Aceasta egalitate arata ca sirul u , u , , um este baza în Im(G) deci Im(G) este spațiu vectorial finit dimensional Astfel, am ajuns la o contradicție si deci G admite un sir infinit de valori proprii încheiem demonstrata ob-servand ca lim Am = +to m^^ Teorema Șirul de funcții proprii {um}m ai operatorului A este un șir ortonormat complet în spațiul Hilbert L (O), iar șirul de funcții proprii um Am este un șir ortonormat complet în spatiul energetic XA m Demonstrația acestei teoreme este laborioasa si nu o facem aici Suntem acum în masura sa formulam urmatoarea teoremă referitoare la solutia generalizată a ecuatiei Au = F, F G L (Q) Teorema Oricare ar fi F G L (Q), solutia generalizată uF a ecuatiei Au = F este dată de: i i uF — T L (Q) ’um , Am unde {um}m este sirul de functii proprii ale operatorului A (A- prelungirea Friedrichs a opertorului A) ortonormal complet în spatiul Hilbert L (Q) Demonstrație: Deoarece uf — l -um sau uf — A= L , (V)v G Xa avem ca: um n||i, = li • l se obtin vectorii bazei ortonormale I nn mn • sin — xi • sin —— x > l J m,n Solutia generalizata este: uF(xi,x ) = m=in=i F (xi, x ) nn mn^ nn mn —j= • sin — xi • sin —— x \ dx-jdxr •sin -— xi-sm-— x vlil li l J ѴІ і li l unde F(xi,x ) este functia data la a), b), c) CAPITOLUL Problema Cauchy-Dirichlet pentru ecuatii parabolice Fie Q C IRn un domeniu marginit cu frontiera dQ neteda (parțial neteda) și functiile reale: aij, c,u : Q IR , f: [ , +x) x Q IR , g : [ , +x) x dQ IR cu urmatoarele proprietati: i) aij sunt functii de clasa C pe Q și aij = aj-i , i, j = , n; c este o functie continua pe Q,u este o functie continua pe Q și de clasa C pe Q ii) exista p > astfel încat pentru orice (^, , £n) G IRn sa aibe loc inegalitatea n n n E E aij(X) • ii • ij > jj E '■ X G Q; i= j= i= iii) c(X) > , (V) X G Q; iv) f este functie continua pe [ , x ) x Q si g este o functie continua pe [ , +x) x dQ Definitia Problema care constă în determinarea funcțiilor reale u : [ , +x) x Q IR care au urmatoarele proprietăți: ) u este continua pe [ , +x) x Q, de clasă C pe ( , +x) x Q si pentru orice t G ( , +x) fixat u este de clasă C pe Q ) du dt +c(X) • u(t,X)=f (t,X), (V) t> §i (V)X G Q ( ) ) u(t,X)= g(t,X), (V) (t,X) G [ , +x) x dQ ( - ) ) u( ,X)= u (X), (V)x G Q ( - ) se numeste Problemă Cauchy-Dirichlet pentru ecuatia parabolica Problema Cauchy-Dirichlet pentru ecuatii parabolice Definitia O funcție u care verifica condițiile din definiția precedentă se numește soluție clasică a Problemei Cauchy-Dirichlet pentru ecuația parabolică Propoziția Dacă exisța un domeniu Q' C IRn care include mulțimea Q si o funcție G : [ , +to) x Q' IR de clasa C pe [ , +to) x Q' asțfel încâț G(t, X) = g(t, X) V(t, X) G [ , +to) x dQ, ațunci Problema Cauchy-Dirichleț neomogenă pențru ecuația parabolica, prin schimbarea de funcție necunoscuța v(t,X) = u(t,X) — G(t,X) se reduce la Problema Cauchy-Dirichleț omogena pențru ecuația parabolică: dv d (dv\ dG dt - E dX ( S aij(X) • ■: ) + C(X) • V(t’XM(t’X) - ~dt + i= \j = j / f fi} ХУ \ +E sx S (X) ■ ăd - C(X)' G(t-X)■ ( - ) (V) t)> si (V)X g Q v(t,X) = (V)(t,X) G [ , +to) x dQ ( - ) v( ,X) = uo(X) — G( ,X) (V) X G Q ( - ) Demonstrație: prin verificare Observația Problema Cauchy-Dirichlet pentru ecuația parabolica cu condiții pe frontiera neomogene (prezentata în Def ), se reduce la o Problema Cauchy-Dirichlet pentru ecuatia parabolica cu conditii pe frontiera omogene Datorita acestui fapt, vom considera în continuare Problema Cauchy-Dirichlet pentru ecuatia parabolica de tipul urmator: CAPITOLUL ъ n n / n p \ d dD* X ' (t’X) i= j= (V)t> si (V) X G ( ) u(t, X) = (V) (t, X) G [ , x ) x d ( ) u( ,X)= uo(X) (V) x G ( ) Daca consideram operatorul diferential A definit pe spatiul de functii: D = {w|w : IR , w G C( ) П C ( ) si w|dQ = }, cu formula: Aw = - + c(X) • w atunci Problema Cauchy-Dirichlet ( )-( ) se scrie: ăt+Au=f ( ) u( , X) = u Problema Cauchy-Dirichlet pentru ecuatii parabolice în continuare, vom formula o problema mai generala pentru care vom demonstra o teorema de existenta si unicitate Teorema Dacă funcția u : [ , +to) x Q IR este soluție clasică a Problemei Cauchy-Dirichlet ( )-( ), atunci funcția V definită prin V (t)(X) = u(t,X) are urmatoarele proprietăți: a) V : [ , +to) L (Q) este continua; b) V : ( , +to) L (Q) este de clasa C ; c) V : ( , +to) XA este continua, unde XA reprezinta spațiul energetic al operatorului A d) V( ) = u dV e) + AV = F(t), (V) t> unde F(t)(X) = f (t,X) Demonstrație: a) Fie t G [ , +to) si n > Functia u(t,X) este uniform continua pe multimea compacta [t — n, t + n] x Q (daca t = , atunci [ , n] x Q) Prin urmare, (V)s > , ( ) (e) > astfel ca (V)(t',X'), (t",X") G [to — n,to + n] x Qcu \t' — "\ , de unde V = V Teorema Daca functia F : [ , +to) L (Q) este de clasa C pe [ , +to) atunci problema Cauchy abstractă pentru ecuatia parabolica are o solutie unică Demonstrație: Va trebui sa aratam doar ca problema ( ) are solutie, unicitatea o avem deja in baza teoremei precedente Presupunem ca avem o solutie si, pentru t G [ , +to) o dezvoltam după sistemul ortonormat complet de functii proprii (um)meIy ale operatorului A : oo V(t) — L (Q) 'um m= Șirul corespunzator de valori proprii va fi notat cu L (Q) •Um m= CAPITOLUL oo V — L (Q) 'um Notam: vm(t) — L (n) fm(t) — L (Q) Vm(t) L (Q) și din ecuația %+AV=F ■ și condiția inițiala V ( ) = Vo deducem: + Am • vm — fm, Vm( ) — V°m m — , , , dt Aceste probleme cu date initiale au solutiile date de formula t Ăm(t-s) • fm(s)ds m — , , , vm(t) — vom e -Amt e de unde rezulta ca soluția V(t) a Problemei Cauchy abstracte ( ) verifica egalitatea: -Amt 'W (t / e-Am(t-s) fm(s)ds I • um ( ) Vom arata acum ca, daca V G L (Q) și functia F : [ , +to) L (Q) este de clasa C pe [ , +to), atunci membrul drept al formulei ( ) definește o functie V(t) care este solutie pentru problema Cauchy abstracta, adica V are proprietatile (a), (b), (c) din Definitia ( ) In prima etapa, trebuie demonstrata convergenta seriilor din membrul drept al egalitatii ( ) si examinata ”netezimea” functiilor care sunt sumele acestor serii Deoarece {um}m este un sistem ortonormat complet în L (Q), din convergen-oo ta seriei numerice V | • e- Xmt rezulta convergenta in L (Q) a seriei de Problema Cauchy-Dirichlet pentru ecuatii parabolice functii E vom • e Xmt m= o • um Seria E |vm| m= •e Xmt (V) t > , este majorata de o seria numerica E |^| care este convergenta (V G L (Q)) Rezulta astfel m= oo ca, seria E in m= spatiul L (Q) si suma ei este functie continua de t o / t \ Pentru a arata ca seria E I f e—Xm(t—s) • fm(s)ds\ • um converge în L (Q) m= \ / uniform pe un segment oarecare [ , T], procedam dupa cum urmeaza: con-oc t sideram seria E | J e-Xm(t—s) • fm(s)ds| , si o majoram astfel: m= e-Xm(t—s) fm(s)ds remarcam ca domeniul D(A) poate fi caracterizat astfel: {сю сю v cmum| Xm • Cm dV Verificarea egalitatilor: ——+ AV = F(t) si V( ) = V este imediata dt Exercițiul Gasiti Problema Cauchy abstracta in cazul Problemei Cauchy-Dirichlet du d u dt a dx t > , x G ( , l) u(t, ) = u(t, l) = u( , x) = u (x) si determinati solutia ”tare” CAPITOLUL Răspuns: bk — k G IV*, uk — sin ^V oo u(x, t) — ak • e-( k=i knx Sin —— unde l ak — knx sin—^dx Exercițiul Determinați soluțiile următoarelor Probleme Cauchy-Dirichlet: du d u — — -, dt dx (x,t) G ( , ) x ( , +x) a) k u(x, ) — sin nx, x G [ , ] R: u(x, t) — • e n t • sin nx du d u t - — •—— + e sin x dt dx x G ( , n) x ( , x ) b) u( , t) — u(n, t) — , t> k u(x, ) — sin x • cos x, x G [ , n] R: u(x, t) — t • e t • sin x + e t • sin x f du d u + x + , (x,t) G ( , ) x ( , +x) c) * u( ,t) — t + , u( ,t) — t + , t > u(x, ) — , x G [ , ] R: u(x, t) — t + + x • t Calculul simbolic si numeric pentru ecuatii parabolice Calculul simbolic si numeric al solutiei Problemei Cauchy-Dirichlet pentru ecuații parabolice Calculul simbolic al soluției unei probleme Cauchy-Dirichlet nu poate fi realizat cu functia pdsolve In astfel de cazuri se trece la rezolvarea numerica Pentru exemplificare, vom considera trei Probleme Cauchy-Dirichlet pentru ecuafii parabolice a caror solutie o vom determina numeric folosind functia pdsolve cu sintaxa pentru calcul numeric: pdsolve(PDE or PDE system, conds, type=numeric, other option); in care: PDEorPDEsystem - ecuatia cu derivate partiale pe care dorim sa o rezolvam sau sistemul de ecuatii cu derivate partiale conds type = numeric otheroption - conditiile initiale si conditiile la limita - indica rezolvarea utilizand metode numerice - diferite optiuni (de ex metoda numerica, nr de puncte, etc ) Exemplul : ' du d u Si (x’ ) = Î • Sx (xt) PDE :=diff(u(x,t),t)= / *diff(u(x,t),x,x); PDE := дu (x,t) = / u (x,t) > IBC := {u( ,t)= , u( ,t)= , u(x, )= }; IBC := {u ( , t) = , u ( , t) = , u (x, ) = } > pds := pdsolve(PDE ,IBC ,numeric); CAPITOLUL pdsl := module () local INFO; export plot, plotdd, animate, value, settings; option ‘Copyright (c) by Waterloo Maple Inc AU rights reserved end module > pl := pdsl:-plot(t= ): p := pdsl:-plot(t=l/ ): p := pdsl:-plot(t=l/ ): p := pdsl:-plot(t=l): p := pdsl:-plot(t= ): plots [display]({pl,p ,p ,p ,p }, title=‘Heat profile at t= , , , , ‘); Heat profile at t= , , , , : : o s; ; ; Qj Figura > pdsl:-plot d(t= ,x= ,axes=boxed); Calculul simbolic și numeric pentru ecuații parabolice Figura Exemplul : — = —— ) (ж,£) + e t • șina? dt \ dx J m( , i) = U (тГ, i) = u(x, ) = cos x sin X Heat equation > PDE := diff(u(x,t),t)= *diff(u(x,t),x,x)+(exp(- *t))*sin(x); PDE := (x, t) = (ж, t) + e~ t sin (ж) > IBC := {u( ,t)= ,u(Pi,t)= ,u(x, )= *cos(x)*sin(x)}; IBC := {м ( , i) = , и (тг, i) = , и (ж, ) = cos (ж) sin (ж)} > pds := pdsolve(PDE ,IBC ,numeric); pdsl := module () local INFO; export plot, plotdd, animate, value, settings; option ‘Copyright (c) by Waterloo Maple Inc AU rights reserved end module > p := pds :-plot(t= ): p := pds :-plot(t=l/ ): CAPITOLUL p := pds :-plot(t=l/ ) : p := pds :-plot(t=l): plO := pds :-plot(t= ): plots [display]({p ,p ,p ,p ,pl }, title=‘Heat profile at t= , , , , ‘); Heat profile at t= , , , , Figura > pds :-plot d(t= ,x= ,axes=boxed); Figura Calculul simbolic si numeric pentru ecuatii parabolice Exemplul : ' du d u — (x, t) = —— (x, t) + t • cosx dt dx PDE := diff(u(x,t),t)=diff(u(x,t),x,x)+t*cos(x); PDE := dPu (x, t) = u (x, t) + t cos (x) > IBC := {u( ,t)=t,u(Pi,t)= ,u(x, )=cos( *x)+cos( *x)}; IBC := {u (n, t) = , u ( , t) = t, u (x, ) = cos ( x) + cos ( x)} > pds := pdsolve(PDE ,IBC ,numeric); pdsl := module () local INFO; export plot, plot d, animate, value, settings; option ‘Copyright (c) by Waterloo Maple Inc All rights reserved ‘; end module > q := pds :-plot(t= ): q := pds :-plot(t= / ): q := pds :-plot(t= / ): q := pds :-plot(t= ): q := pds :-plot(t= ): plots[display]({q ,pq ,q ,q ,q }, title=‘Heat profile at t= , , , , ‘); CAPITOLUL > pds :-plot d(t= ,x= ,axes=boxed); Figura Problema Cauchy-Dirichlet pentru ecuatii hiperbolice Problema Cauchy-Dirichlet pentru ecuatii hiperbolice Fie Q C Rn un domeniu marginit cu frontiera dQ neteda (parțial neteda) și funcțiile reale aij, c, u , u : Q IR , f : [ , +to) x Q IR , g : [ , +to) x dQ IR cu urmatoarele proprietati: i) aij sunt functii de clasa C pe Q și aij = aj-i, i, j = , n; c este continua pe Q; u este de clasa C pe Q și de clasa C pe Q; u este continua pe Q si de clasa C pe Q ii) exista p > astfel încat pentru orice £n) G RRn sa aibe loc inegalitatea: n n n aij (X) • ^^ & (V)x G Q i= j= i= iii) c(X) > , (V)X G Q; iv) functiile f : [ , +to) x Q Ri si g : [ , +to) x dQ Ri sunt continue Definitia Problema care consta în determinarea funcțiilor reale u : [ , +to) x Q R continue pe [ , +to) x Q, de clasa C pe [ , +to) x Q si de clasă C pe ( , +to) x Q, care au următoarele proprietăți: r\ П r\ f П O \ d u v—л d /v—л , x du \ /-r^\ p / t^x да - L ■ aij (X) • ăd+c(X) •u=f (t’X ’■ i= j= ( ) (V)(t, X) G ( ,+to)x Q u(t,X) = g(t,X), (V)(t,X) G [ , +to) x dQ ( ) u( ,X)= uo(X), (V)X G Q ( ) du — ( ,X)=u (X), (V)X G Q ( ) se numește Problemă Cauchy-Dirichlet pentru ecuatia hiperbolică CAPITOLUL Definitia O funcție u care verifica condițiile din definiția precedentă se numește soluție clasică a Problemei Cauchy-Dirichlet pentru ecuația hiperbolică Propoziția Dacă exisța un domeniu O' C IIC care include domeniul O si o funcție G : [ , +to) x O' IR de clasă C pe [ , +to) x O' asțfel încâț G(t, X) = g(t, X), (V)(t, X) G ( , +to) x dO ațunci Problema Cauchy-Dirichleț pențru ecuația hiperbolică, prin schimbarea de funcție necunoscuță v(t,X) = u(t,X) — G(t,X), se reduce la Problema Cauchy-Dirichleț pențru ecuația hiperbolicăa cu condiții nule pe fron-țieră: n d v dt д I dv \ • dr) + c(X)’v(t,X) = г \j=i j/ ( ) = f (t,X) d G dt — c(X) • G(t,X), (V)(t, X) G ( , +to) O x v(t, X) = , (V)(t,X) G [ , +to) x O ( ) u( ,X) = uo(X) — G( ,X) = vo(X), (V)X G O ( ) ( ) dv dG — ( ,X) = ui(X) — — ( ,X) = Vi(X), (V)X G O Demonstrație: Prin verificare Observația Propozitia reduce problema ( - ) cu conditie nenula pe frontiera la problema ( - ), in care conditia la frontiera ( ) este zero: v(t, X) = , (V) (t, X) G [ , +to) x dO De aceea, vom studia existența și unicitatea soluției Problemei Cauchy-Dirichlet pentru ecuatia hiperbolica cu conditia la frontiera nula Problema Cauchy-Dirichlet pentru ecuatii hiperbolice Vom considera în continuare Probleme Cauchy-Dirichlet pentru ecuația hiperbolica de urmatoarea forma: d u \—л d I \—л , x du \ , x , x , x ătr-EădEaii(X)• ăd+ ' >• u(t-X '(t-X>■ i= г \j= -v ( ) (t, X) G (Q+to)x Q u(t, X) = , (V)(t, X) G [ , +to) x Q ( ) u( , X) = uo(X), (V)x G Q ( ) dt( ,X)= ui(X), (V)X G Q ( ) în care functiile aij, c, f au proprietatile (i), (ii), (iii), (iv) anterior prezentate; functia u este de clasa C pe Q si de clasa C în Q; functia u este continua pe Q si de clasa C în Q O solutie clasica a acestei probleme este o functie u : [ , +to) x Q IR care este de clasa C pe [ , +to) x Q, si este de clasa C pe ( , +to) x Q si verifica: rj n fi) ( n o \ X (X)• ău) + c(X)• u(t X)=f (t,X); (t,X)G i= i j= j u(t, X) = , (V) (t, X) G [ , +to) x Q u( , X) = u (X), (V)X G Q dt( ,X) = U (X), (V)X G Q Considerând operatorul diferential A definit pe spatiul de functii D = {w|w : Q IR ; w G C(Q) C C (Q), W|dQ = } CAPITOLUL cu formula Aw = ecuația hiperbolica poate fi scrisa sub forma: — + A • u(t, X) = f (t, X), (V)(t, X) G ( , +to) x Q dt Teorema Dacă funcția u : [ , +to) x Q IR este soluție clasică a problemei ( - ), atunci funcția U definită prin: U (t)(X) = u(t,X) are urmatoarele proprietăți: i) U : [ , +to) L (Q) este de clasa C pe [ , +to) ii) U : [ , +to) H) este functie continua; iii) U( ) = u si U'( ) = u Demonstrație: i) Aratam la început ca functia U : [ , +to] L (Q) este derivabila în orice t G [ , +to) Pentru aceasta, consideram t G [ , +to) si apoi raportul -ă-[U(t) - U(to)] G L (Q) t — t si aratam ca acest raport tinde la — (t ,x) G L (Q) în norma L (Q) atunci cand t t Aceasta înseamna ca trebuie sa arătam egalitatea: lim [ T - [U(t) — U(С)] (X) Q t — t dX = dt ' ) Folosind teorema creșterilor finite a lui Lagrange scriem: — — [U(t) — U(t )] (X) = - — [u(t, X) — u(t , X)] t — t t — t =ă(t(X ’-X) Problema Cauchy-Dirichlet pentru ecuatii hiperbolice cu \t(X) - to| ) si prin urmare: (V)s > , ( ) (s) a î (V)(tz, X'), (t”, X”) G [t — n, t + n] x Q cu [t' — t” | * Pentru aceasta, vom folosi egalitatea dt ') IIU'(t) — U'(to)Hb («) = [ Jq du du ^(t,X) (to ,X) dX dt dt CAPITOLUL și faptul ca funcția — este continua pe [ , +to) x Q Din continuitatea dt du functiei — rezulta ca aceasta este uniform continua pe o multime de forma dt [t — n,t + n] x Q, (y > ) și prin urmare: (V)e > , ( ) (s) a î (V)(t/, X'), (t”, X ”) G [t — n, t + n] x Q daca |t — t”| ’x > |t — t | dX+ d u BtBxi (t*(X ),X) |t — t | dX Problema Cauchy-Dirichlet pentru ecuatii hiperbolice CU \t(X) — t | astfel ca \\u(t) — U(to)HHox si subspatiul de functii ST definit prin: St = {V G S\V(T) = } Teorema Dacă funcția u = u(t,X) este o soluție clasică a problemei ( - ), atunci funcția U definita prin U(t)(X) = u(t,X) aparține spațiului de funcții S, verifica U( ) = u , U'( ) = u si pentru orice T > si orice funcție V G ST are loc egalițațea: T У (t),V>(t) >L (Q) dt— L (q) + ( ) T T У A dt L (Q) dt Demonstrație: Apartenenta functiei U la spatiul S a fost demonstrată S-a arătat de asemenea că U( ) = u si U'( ) = u Rămăne doar să arătăm că CAPITOLUL pentru orice T > si orice V G ST are loc egalitatea ( ) In acest scop, fie T > si V G ST Scriind egalitatea ( ) sub forma: ^^(X) + A • U(t)(X) = F(t)(X) (unde F(t)(X) = f (t,X)), înmulțind aceasta egalitate cu funcția V(t)(X) și integrând pe Q obținem: d U L W + a= L (Q) • Integram acum aceasta egalitate în raport cu t pe segmentul [ , T], tinem seama de V (T) = si obtinem: L (n) |T t T -j l (q) dt+У a dt = oo = [ l (q) dt Jo adica: T T L (Q) L (Q) dt+У A dt = oo = [ L (Q) dt o ceea ce este echivalent cu: t T l (Q) — У L (Q) dt + У A= oo = [ l (q) dt o Problema Cauchy-Dirichlet pentru ecuatii hiperbolice Definiția O funcție U E S se numește soluție generalizată a problemei ( - ) dacă U( ) = u , U'( ) = u și pentru (V)T > , (V)V E ST verifică: T — l (q) — I l (q) dt+ + [ A dt = [ L (Q) dt ( - ) Observația O soluție clasica u(t,X) a problemei ( - ) definește o soluție generalizata a acestei probleme Teorema Dacă funcția U E S este soluție generalizată a problemei ( - ) si functia u : [ , +to) x Q U’ definită prin u(t,X) = U(t,X) este de clasa C pentru t > si X E Q, atunci functia u = u(t,X) este solutie clasicăa a problemei ( - ) Demonstrație: Egalitatea ( ): u(t,X) = , (V)t > si (V)x E dQ rezulta din apartenenta U(t) E H Egalitatea ( ): u( ,X) = u ( ), (V)X E Q rezulta din U( ) = u du Egalitatea ( ): — ( ,X) = u (X), (V)X E dQ rezulta din U'( ) = u Ramane doar sa aratam egalitatea ( ) adica: — + A • u(t,X) = f (t,X) Pentru a deduce aceasta egalitate pornim de la egalitatea ( ) pe care o scriem sub forma: T — [ ( [ (t,X) • V'(t)(X)dX]dt — [ ui(X) • V( )(X)dX+ J \ Jq dt J Jq CAPITOLUL A • u(t,X) • V(t)(X)dX T f (t,X) • V(t)(X)dX dt Itervertind ordinea de integrare și făcând o integrare prin părți in raport cu t în prima integrala, obtinem: T — + A • u(t,X) - f (t,X) • V(t)(X)dtdX = pentru orice V G ST (am tinut seama de faptul ca V(T) = ) Rezulta de aici ca: — + A • u(t,X) = f (t,X), (V)t> , (V)X G Q dt Observația Aceasta teorema arata ca dacă o solutie generalizată este suficient de ”neteda” atunci ea este solutie clasica Teorema (de unicitate a soluției generalizate) Problema (l )-(‘ ) are cel mult o soluție generalizată Demonstrație: Presupunem prin absurd ca problema ( )-( ) admite solutiile generalizate U (t) si U (t) si consideram functia U(t) = U (t) — U (t) Pentru orice T > si orice V G ST functiile U si V verifica: T T У L (^) dt + j a dt = T Functia V definita prin V(t) = f U(t)dr apartine spatiilor de functii S si ST t si are urmăatoarele proprietăati: V'(t) = —U(t) si V"(t) = — U'(t) înlocuind acestea în relatia de mai sus rezulta că: T T J l^) dt — f A dt = Problema Cauchy-Dirichlet pentru ecuatii hiperbolice Deoarece — — d — L (q) = • — || V (t) ||L (Q) si — — d — A= - • T|| V(t)||A A dt" v egalitatea precedenta implica egalitatea: - HV'( )||L (n) - \\v(T)||A + IIV( )||A = "Ținem seama acum de egalitățile V(T) = , V'( ) = U( ) = și deducem ca: \\V'(T )ПІ (П) + \\V'(t )||A = , din care rezulta V(T) = și V( ) = Intrucat T > este oarecare rezulta V(T) = —U(T) = Prin urmare Ui(T) = U (T), (V)T > Astfel, rezulta cele doua solutii generalizate coincid Consecința Problema ( )-( ) are cel mult o soluție clasică Teorema (de existență a soluției generalizate) Dacă funcția F definita prin F(T)(X) = f (t,X) este continuă ca funcție cu valori în L (Q) și dacă u G H , ui G L (Q), atunci problema ( )-( ) are o soluție generalizată Demonstrație: Facem demonstrata în doua etape In prima etapa deducem o formula de reprezentare a solutiei generalizate in ipoteza ca aceasta solutie exista In a doua etapa aratam ca formula de reprezentare gasita în prima etapa, în conditiile teoremei, defineste o functie care este o solutie generalizata Etapa I Presupunem ca U = U(t) este o solutie generalizata a problemei ( )-( ) si consideram sirul valorilor proprii L (Q) CAPITOLUL și reprezentarea OO U(t) У Uk(t) • Uk k= a funcției U(t) Pentru ca U G C ([ , +to); L (Q)) funcțiile uk(t) sunt deriv-abile si au derivata continua, iar U'(t) se reprezinta astfel: OO U/(t) = J Uk(t) Uk• k= In virtutea acestei formule de reprezentare, egalitatea ( ) pe care o satisface solutia generalizata U devine: T ^ Uk(t) L (Q) dt — uk( )^ L (Q) + k= k= TT p Ak • Uk(t) L (Q) dt = Y fk(t) L (Q) dt k= k= Fie acum j un numar natural oarecare fixat si functia Vj G ST definita prin: Vj(t) = (T — t) • Uj • In egalitatea precedenta înlocuim V cu Vj, tinem seama de egalitatile l (q)= dkj; V'(t) = — Uj, V( ) = T • Uj si obtinem: T T T —T• l (q) + j uj(t)dt + Aj • !(T — t) • fj(t)dt = f (T — t)fj(t)dt pentru j = , , • • • Derivand de doua ori in raport cu T rezulta: u'j (T) + Aj • Uj (T) = fj (T), (V)T> uj( ) = L (Q) Uj( ) = L (Q), j = , , • • • ( ) Problema Cauchy-Dirichlet pentru ecuatii hiperbolice Problema cu datele inițiale ( ) are o singură soluție și aceasta este data de: u-(t) — , j — , , o Astfel rezulta ca solutia generalizată U are urmatoarea reprezentare: и (t) —E j= L (Q) • cos y/Xj • t • £ (Qysin y/Xj • t • Uj + +E j= T h (T ) o • sin Xj • (t — • Uj + ( ) Etapa II Aratam acum ca, în conditiile din teorema formula ( ) de-fineste o functie U care apartine spatiului S si verifica ( ) pentru T > Convergentele У^ | L | L | £ (q) • cos у/Xj • t +-• l • sin y/Xj • t Xj Uj si faptul ca suma seriei este functie continua de t cu valori în L (Q) Inegalitatile: T У fj(t) • sin y/Xj • (t — t)dT o T “ Ti I fj (T)dT, o CAPITOLUL (V)t G [ ,T ],j = , , precum și convergența seriei de funcții continue și pozitive У^ ff(r) la funcția j=i conținua \\F(T)|| (n) implica convergența uniforma în raport cu t G [ , T] ((V)T > și T £ (q) • sin • t+ £ (Q) • coș у/Xj • t • Wj+ j= Wj T o și se arața ca aceasța converge uniform în raporț cu t G [ ,T](T > și T l (q) | esțe convergența Aceasța din urma, poațe fi j=i Problema Cauchy-Dirichlet pentru ecuatii hiperbolice scrisa sub forma: У^ Aj • I l (q) | У V ’ ’ I l (Q) | — j=i j=i Aj j=i si prin urmare este convergenta pentru ca u G H Următoarea serie a cărei convergenta trebuie examinata este seria £ (q) • cos у/Aj • t • Uj j=i Aceasta serie este uniform convergenta în raport cu t G [ , +to) daca seria numerica | L (q ) | este convergenta , ceea ce este adevarat j=i pentru ca u G L (Q) Ceea de a treia serie care trebuie examinata este seria: fj (t) • cos y/Xj • (t — t)drUj Aceasta este uniform convergenta in raport cu t G ( , T] (T > si T L (Q) • cos \ Ă • t • Uj j= în H poate fi asigurata prin convergența uniformă în raport cu t G [ , +to) a seriei: У VĂ" L (Q) • cos • t • —= în H j=i VAj Convergența acestei serii rezulta din convergența seriei numerice У? ĂjI L (Q) | , j=i convergența ce este asigurata de ipoteza u G H Convergența uniforma în raport cu t G [ , +to) a seriei de funcții i i У^ —Д= • L (Q) • si^/Ăj • t • Uj în spațiul H este asigurata de convergența seriei У^ | |L (q), j=i convergența care rezulta din apartenența u G L (Q) Problema Cauchy-Dirichlet pentru ecuatii hiperbolice In fine, convergența uniforma în raport cu t G [ , T] (T > și T n u(x, ) = sin —x, [ ,l] â u â u ât âx R: u(x, t) t • sin x, [ ,l] n n = cos l t • sin l x (x, t) G ( ,n) x ( , +to) u( , t) = u(n, t) = , t> u(x, ) = sinx cos x, x G [ , n] Й ■ x G x G — (x, ) = sin x, О/ ’ R: u(x, t) d u d u = , o/ dx x G [ , n] = (sin t — t) • sin x + cos t • sin x+ +- sin t • sin x (x,t) G ( ,n) x ( , +to) u( , t) = u(n, t) = , t> u(x, ) = sinx, x G [ , n] — (x, ) = sin x + sin x, x G [ , n] R: u(x, t) = ( cos t +— sin t | • sin x + - sin t • sin x Calculul simbolic si numeric pentru ecuatii hiperbolice Calculul simbolic si numeric al solutiei Problemei Cauchy-Dirichlet pentru ecuații hiperbolice Deoarece pentru calculul simbolic al soluției unei probleme Cauchy-Dirichlet funcția pdsolve nu afiseaza nimic, vom trece la rezolvarea numerica folosind functia pdsolve specifica calculului numeric, a carei sintaxa a fost prezentata într-unul din paragrafele anterioare Pentru exemplificare consideram urmatoarele probleme Cauchy-Dirichlet de tip hiperbolic: Exemplul : d u d u аё (x,t) = (x,t) u( , t) = u( , t) = PDE :=diff(u(x,t),t,t)=diff(u(x,t),x,x); PDE := u (x,t) = u (x,t) > IBC := {u( ,t)= , u( ,t)= ,u(x, )=x~ -x, D[ ](u)(x, )= }; IBC := {u ( , t) = , u ( , t) = , u (x, ) = x — x, D (u) (x, ) = } > pdsl := pdsolve(PDE ,IBC ,numeric); pdsl := module () local INFO; export plot, plot d, animate, value, settings; option ‘Copyright (c) by Waterloo Maple Inc All rights reserved ‘; end module > p := pds :-plot(t= ): CAPITOLUL p := pdsl:-plot(t=l/ ): p := pdsl:-plot(t=l/ ): p := pdsl:-plot(t=l): p := pdsl:-plot(t= ): plots[display]({pl,p ,p ,p ,p }, title='Wave profile at t= , , , , ‘); Figura > pdsl:-plot d(t= ,x= ,axes=boxed); Figura Calculul simbolic si numeric pentru ecuatii hiperbolice Exemplul : d u d u ^ (x,t) = dx (x,t) - t Sinx u( , t) = u(n, t) = u(x, ) = • sinx • cosx — (x, ) = • sin x Wave equation > PDE :=diff(u(x,t),t,t)=diff(u(x,t),x,x)-t*sin(x); PDE := d u (x, t) = -^ u (x, t) — t sin (x) > IBC :={u( ,t)= ,u(Pi,t)= ,u(x, )= *(sin(x))*(cos(x)), D[ ](u)(x, )= *sin( *x)}; IBC := {u( , t) = , u(n, t) = , u(x, ) = sin(x) cos(x), D (u)(x, ) = sin( x)} > pds := pdsolve(PDE ,IBC ,numeric); pdsl := module () local INFO; export plot, plot d, animate, value, settings; option ‘Copyright (c) by Waterloo Maple Inc All rights reserved ‘; end module > p := pds :-plot(t= ): p := pds :-plot(t= / ): p := pds :-plot(t= / ): p := pds :-plot(t= ): p := pds :-plot(t= ): plots[display]({p ,p ,p ,p ,p }, title=‘Wave profile at t= , , , , ‘); CAPITOLUL Figura > pds :-plot d(t= ,x= Pi/ ,axes=boxed); Figura Calculul simbolic si numeric pentru ecuatii hiperbolice Exemplul : d v d v ăl PDE :=diff(v(x,t),t,t)= *diff(v(x,t),x,x); PDE := v (x,t)= djrv (x,t) > IBC :={v( ,t)= ,v(Pi/ ,t)= ,v(x, )= , D[ ](v)(x, )= *sin(x); IBC := {v ( , t) = , v ( / n, t) = , v (x, ) = , D (v) (x, ) = sin (x)} > pds := pdsolve(PDE ,IBC ,numeric); pdsl := module () local INFO; export plot, plot d, animate, value, settings; option ‘Copyright (c) by Waterloo Maple Inc All rights reserved ‘; end module > q := pds :-plot(t= ): q := pds :-plot(t= / ): q := pds :-plot(t= / ): q := pds :-plot(t= ): q := pds :-plot(t= ): plots[display]({q ,pq ,q ,q ,q }, title=‘Wave profile at t= , , , , ‘); CAPITOLUL Figura > pds :-plot d(t= ,x= Pi/ ,axes=boxed); Figura Bibliografie [ ] V I Arnold, Ecuații diferențiale ordinare, Editura Științifică și Enciclopedica, București, [ ] Șt Balint, A M Balint, S Birăuaș, C Chilărescu, Ecuații diferențiale și ecuații integrale, Editura Universitari de Vest, [ ] V Barbu, Ecuații diferențiale, Editura Junimea, Iași, [ ] R Cristescu, Elemențe de analiză funcțională și ințroducere în țeoria disțrubuțiilor, Editura Tehnica, [ ] S A Coddington, N Levinson, Theory of ordinary differențial equațions, McGraw Hill, New York, [ ] A Corduneanu, Ecuații diferențiale cu aplicații în elecțroțehnică, Editura Facla, , Timisoara [ ] C Corduneanu, Ecuații diferențiale și ințegrale, Universitatea ”Al I Cuza”, Iasi, [ ] G M Fihtenholt, Curs de calcul diferențial și ințegral, Editura Tehnica, [ ] D Gaspar, N Suciu, Analiza complexă, Editura Academiei Romane, [ ] A Halanay, Ecuații diferențiale și ințegrale, Universitatea din Bucuresti, [ ] A Halanay, Ecuații diferențiale, Editura Didactica si Pedagogica, Bu-curesti, BIBLIOGRAFIE [ ] J Ray Hanna, John H Rowland, Fourier series, transforme, and bound-ary value problems A Wiley-Interscience Publication John Wiley and Sons, Inc , , USA [ ] D Haragus, Ecuații cu derivate parțiale, Editura Universitații de Vest, Timișoara, [ ] A Eckstein, D Haragus, Exerciții standard de ecuații cu derivate parțiale, Tipografia Universitari de Vest din Timisoara, [ ] Ph Hartman, Ordinary Differential Equations, J Wiley, New York, [ ] M Hirsh, S Smale, Differential Equations, dynamical systems and linear algebra, Academic Press, New York, [ ] J Hubbard, B West, Equations differentielles et systemes dynamiques, Cassini, Paris, [ ] St Mirica, Ecuații diferențiale, Editura Universități din Bucuresti, [ ] Gh Morosanu, Ecuații diferențiale Aplicații, Editura Academiei, Bucuresti, [ ] M Reghis, P Topuzu, Ecuații diferențiale ordinare, Editura Mirton, [ ] E Rogai, Exerciții și probleme de ecuații diferențiale și integrale, Editura Tehnica, [ ] N Rouche, P Habets, M Laloy, Stability Theory by Liapunov direct Method, Springer, [ ] I A Rus, P Pavel, Ecuații diferențiale, Editura Didactica si Pedagogica, Bucuresti, [ ] I A Rus, Ecuații diferențiale, ecuații integrale si sisteme dinamice Casa de Editura Transilvania Press, [ ] V V Stepanov, Curs de ecuații diferențiale, Editura Tehnica, Bucuresti, BIBLIOGRAFIE [ ] M Stoka, Culegere de probleme de funcții complexe, Editura Tehnica, [ ] M Tihonov, A A Samarski, Ecuațiile fizicii matematice, Editura Tehnica, [ ] K Yosida, Equations differentielles et integrales, Dunod, Paris, https://neculaifantanaru com/calitatile-unui-lider html https://neculaifantanaru com/en/qualities-of-a-leader html